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a grand-pere 



Abstract 



To examine the loss of entanglement in a two-particle Gaussian sys- 
tem, we couple it to an environment and use the Non-Rotating Wave master 
equation to study the system's dynamics. We also present a derivation of 
this equation. 

We consider two different types of evolution. Under free evolution we find 
that entanglement is lost quickly between the particles. When a harmonic 
potential is added between the particles, two very different behaviours can 
be observed, namely in the over and under-damped cases respectively, where 
the strength of the damping is determined by how large the coupling to the 
bath is with respect to the frequency of the potential. 

In the over-damped case, we find that the entanglement vanishes at even 
shorter times than it does in the free evolution. In the (very) under-damped 
case, we observe that the entanglement does not vanish. Instead it oscillates 
towards a stable value. 



Introduction 



Freud said that mankind's ego underwent three downfalls caused by 
Galileo's heliocentrism, Darwin's theory of evolution and his own theory of 
the unconscious mind. Similarly, it can be said that the field of physics 
was "reduced" three times with Newtonian mechanics, statistical mechan- 
ics (followed by quantum mechanics) and the theory of relativity. Whereas 
Newtonian mechanics concerns solid objects, statistical mechanics concerns 
ensembles of particles and is therefore much better suited to study the very 
small. Finally quantum mechanics allows for a study of particles and the 
evolution of their states in time. Entanglement is a property very specific 
to quantum mechanics and is proving a wonderful resource to the field of 
quantum information and quantum computing. This work sets to study how 
a quantum state, and the entanglement within it, evolve with time. 



1.1 Study of open quantum systems 

The mechanics of closed quantum mechanical systems are rather well- 
known and have been extensively studied. If one considers a system with 
wavefunction \I/(r, t), the time-evolution of the system is obtained from the 
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Schrodinger equation 



H\V(r,t))=ih-\*(r,t)) 



(1.1) 



where H is the Hamiltonian of the system, usually in terms of the Lapla- 
cian V, a potential V(r), the mass of the particle m and Planck's constant. 



When the particle is coupled to an environment (which in the remainder 
of this work will also be known as a heat-bath or a reservoir), the Schrodinger 
equation is no longer sufficient to describe the system's evolution, since it fails 
to take the bath's variables and behaviour into account. The Schrodinger 
equation described above has to be modified to account for the damping 
caused by the environment. It is now appropriate to examine the relationship 
between the bath and the system and the dynamics of the dissipation from 
a general point of view. We now consider one-dimensional systems. 

If the system is interacting with the bath via a random force F(t), 
the system evolution is called stochastic. Classically, a stochastic equation 
may be written for a physical variable q and the random force is real- valued. 
Quantum mechanically, the variable q becomes an operator and the random 
force an operator-valued function. One can write a stochastic equation as 



H = -— V 2 + V(r) 



(1.2) 



q + 1 q + uj 2 q = F(t)/m (1.3) 
A more general form of this equation is the Quantum Langevin Equation 
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(Q.L.E.) |L088 



mx + I dt'n(t - t')x(t') + V'(x) = F(t) (1.4) 

J — oo 

where the dot and the prime are respectively the derivative with respect 
to t and x. The function fi(t) is a "memory" function and depends on the 
way the bath is coupled to the system. It can be reduced to a constant 7 in 
a Markov approximation. The QLE can be derived from the Hamiltonian of 
the complete system using the Heisenberg equations of motion and a specific 
bath model (see Chapter 3). 

The state of a quantum system is best described by a density matrix p. 
A state is called pure when its density matrix is the projection | \1/ ) ( \1/ | onto 
a vector in the Hilbert space. The time evolution of the complete system 
(system + bath) may then be expressed by the Von Neumann equation 



with a Hamiltonian chosen depending on the system and the type of cou- 
pling. Considering only the density matrix of the system, one may derive a 



macroscopic description in terms of the operators x and p from (1.5), which 
is generally known as a master equation. There are various master equations 
in common use. The types of master equations depend on the assumptions 
made concerning the bath, its coupling to the system and also the initial 



conditions. A general derivation may be found in LO96,GZ00 . 

A very popular type is obtained in the rotating-wave approximation 
(RWA). The main reason for the popularity of this approximation resides 



in the fact that a master equation obtained in the RWA is of the following 
Lindblad form Lin 76 



where H is the Hamiltonian and Aj are operators depending on the system 
studied. The solutions of Lindblad type equations are guaranteed to be 
positive semi-definite matrices. The RWA is well suited for systems with 
well-defined energy levels and rather weak coupling, for instance a harmonic 
oscillator in a weak radiation field. 

One may wish to study more general systems. The coupling may be 
stronger, or the system may have energy levels much closer together. The 
approximations required in this case are quite different and one may de- 
rive a "pre-Lindblad" type of master equation. Various such equations have 
been derived using a variety of approaches. For instance, the Hu-Paz-Zhang 



equation was derived using path integral techniques PZ92], whereas Agarwal 



derived a general equation in phase space formalism Aga71, Aga69, AW68 



An alternative derivation of the Hu-Paz-Zhang equation was proposed by 



Halliwell and Yu HY96| using the Wigner function, a preferred method of 



study for most of these authors, mainly because it is real everywhere. The 



Hu-Paz-Zhang was also solved explicitly by Ford and O'Connell FOOlb and 
the solution was shown to be valid only for high temperatures. The solution 
of a pre-Lindblad equation is not always positive semi-definite and one must 
check in each case whether the resulting solution is a density matrix. 

The present work will use the so called Non- Rotating Wave (NRW) master 
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equation IGZ00LIMG96 



% \H n(+\\ 



2h 



[X, [X,p(t)}. 



[X,[X,p(t)]} (1.7) 



with X = l[H syst , X] and X in principle an arbitrary system operator coupled 
to the bath. It is of pre-Lindblad nature and obtained in the high temper- 
ature limit. It is similar in its formulation to that derived by Agarwal and 



used by Savage and Walls in SW85b,SW85a . The fourth chapter of this 



work contains an analysis of this equation for a single particle with results 
similar to those of Savage and Walls. This equation can be derived from 
the Quantum Langevin Equation for an independent-oscillator model of the 
heat-bath. In the third chapter of this work, such a derivation is presented 
using perturbation theory techniques. 



1.2 Entanglement 

Entanglement is one of quantum mechanics' most fascinating features. 
It is the property two quantum systems can share that allows one to get 
a piece of information about both systems while measuring it in only one. 
It was highlighted by the Einstein, Podolosky and Rosen (EPR) paradox 



in their famous 1935 paper PR35 . The authors used a thought experi- 
ment to express their doubt that the wavefunction description is complete. 
Performing a measurement reduces the system's wavefunction. Performing 
a measurement on one of two correlated systems reduces the two systems' 
wavefunctions simultaneously. According to special relativity, information 
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cannot travel faster than light, yet the systems can be infinitely far apart. 
Such correlation is known as entanglement, a term originally introduced by 



Schrodinger Sch35 



By its very nature, entanglement has proven to be a fantastic resource for 
quantum information. It is used extensively in quantum information theory, 
quantum coding and quantum cryptography, because it allows for very dense 
coding and totally secure encryption. 

Decoherence versus Entanglement Sudden-Death 

When subjected to an environment, a quantum state tends to decohere. 
Decoherence can be qualitatively defined as the destruction of the interfer- 



ence pattern of a quantum state FOOla, O'C05 FO03, LO01, FO04 . This 



definition can be quantitatively represented in the interference term of the 
probability distribution of a particle coupled to an arbitrary reservoir for 
instance. Furthermore, O'Connell O'C05| points out that one may obtain 



different results depending on various conditions, for instance whether one 
assumes initial coupling between the particle and the reservoir or if external 
forces come into play. The decoherence time is typically much shorter than 
the relaxation time of the overall system and as such is a rather important 
quantity Zur91 Zur03| . The degree of decoherence is represented by the 



variance of the off-diagonal terms of the density matrix. 

If one studies the evolution of the entanglement in a two-particle sys- 
tem coupled to two independent reservoirs, one finds that the entanglement 
measures typically show a sharp decrease, vanishing at a finite time whereas 
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the coherence merely vanishes asymptotically. Such phenomenon is known as 
entanglement sudden-death and has been the object of much study in recent 



years |YE03[|YE04[|YE06[|PE01|RM06| . One can typically observe this decay 
by evolving the initial state with a master equation then estimating the de- 
gree of entanglement using an entanglement measure. Eberly et al. uses the 
concurrence 



KB05,Woo98| whereas Eisert et al. observed entanglement de- 



cay and entanglement transfer with the logarithmic negativity HE04J. The 
following subsection explores some entanglement measures in more detail. 
Entanglement arising between the two reservoirs (known as sudden birth of 



entanglement) has also been observed FT06,FT08 ,SR08 . One finds that it 



increases as the entanglement between the two particles decreases. 



Entanglement Measure 

Entanglement of a general quantum state is defined as the opposite of 
separable. Separability can be expressed simply as follows 

Definition 1 A state p of a composite system of two parts is separable if it 
can be written as 

n 

Pab = Y,Pa®Pb ( L8 ) 

k=0 

When the state under study is that of a pure bipartite state Pab, the Von 
Neumann entropy of the reduced density operator, also called the entropy of 
entanglement, is a good measure of entanglement. It is the quantum analogue 
of the Shannon entropy and is defined by S(pa) = —T^(pa ^(pa)), where 
Pa = Tr^(p). The Von Neumann entropy is invariant under a closed system 
time evolution, as will be proven later in this work. 
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When studying mixed systems or systems with more than two particles, 
the Von Neumann entropy is no longer sufficient to fully quantify the entan- 
glement. To be called an entanglement measure, a function £ must satisfy 



the following three conditions RK97 



1. E(p) vanishes if the state p is separable. 

2. £ does not increase on average under local operations and classical 
communications (LOCC). 

3. £ is invariant under local unitary transformations. 

These conditions however do not uniquely specify a measure for mixed 
states. Numerous measures have been proposed, such as the entanglement of 
formation Woo98 , the entanglement of distillation SW96|, the concurrence 



KB05, Woo98 or the global entanglement for multipartite system MW02 



Ved08 , to name but a few. The concurrence in particular has been quite 



popular YE03,YE04 ( YE06 PE01 RM06 



One may write the entangled state pab as p = YliPi II II- The en- 

tanglement of formation is the amount of entanglement needed to create the 
entangled state p and is defined as £(p) = min YliPi^(Pi)- The entanglement 
of distillation is the amount of entanglement that one obtains after purify- 
ing the state. The concurrence is related to the entanglement of formation 
and provides a formula for the an abritrary state of two-qubits. It is given 
as C(p) = max{0, Ai — A2 — A3 — A4} where the Aj are the eigenvalues, in 
decreasing order of the Hermitian matrix R = \J yfppyfp- The global entan- 
glement for multipartite systems provides the entanglement of one of the par- 
ticles to the rest of the system and can be written as £ (p) = 2 — -j| Ylj=i Trp^. 
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Chosen here for the ease of computation with which it is computed for 
continuous variables and its widespread use is the negativity, seconded by its 



close cousin, the logarithmic negativity VW02 . They are based on the trace 



norm of the partial transpose of the density operator p where p represents 
a generic state of a bipartite system. The partial transpose is obtained by 
VW02[|And03j|PE03[|HE04l 



(l A 



3b I P Ta I k A ,l B ) = (k A ,j B | p | i A ,l B ) 



The trace norm of a Hermitian operator A is || A \\i= TrV A^A = |A$| 
where the A, are the eigenvalues of A. Density matrices are Hermitian matri- 
ces and as such have positive eigenvalues and || p \\i= Trp = 1. The partial 
transpose p TA also has trace 1 but since it may have negative eigenvalues p^, 
its trace norm reads 



|| P Ta || 1= 1 + 2 | \= l + 2N{p) (1.9) 

i 

where N(p) is the negativity. One can write the partial transpose sepa- 



rability criterion as VW02.KB05 



Theorem 1 Let p TA be the partial transpose of a state p with eigenvalues \. 
If one of the \ is negative, then the state is entangled. 

However, the converse has been shown to be true only for 2x2 and 2x3 
dimensional systems KB05 , S"im00| , i.e. in systems of this size, the lack of 
a negative eigenvalue is not enough to guarantee that the state is separable. 
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The negativity of a separable state p s can be shown to be M{p s ) 
also monotonous under LOCC. 

The logarithmic negativity is expressed as 



= 0. It is 



CM=\og 2 \\fF*\\ x (1.10) 
and has an interpretation as an asymptotic entanglement cost, which 



itself is the asymptotic version of the entanglement of formation PE03 



Since N(p s ) = for a separable state, it is easy to see that Ctf(p) = also. 



Motivation 

Entanglement is a remarkable resource. It is, however, rather fragile 
and examining just how fragile is the aim of this work. To this end, the 
dynamics of the entanglement in a bipartite Gaussian state will be examined. 

The state is prepared such that it is initially entangled. Using the Von 
Neumann entropy, the second chapter will establish that in a closed system, 
this entanglement does not vary. The second part of this chapter will present 
the covariance matrix formalism that will be used in subsequent chapters. 

The state is then subjected to an open system evolution and the Non- 
Rotating Wave master equation is solved using two Hamiltonians. The first 
Hamiltonian considered is a free-particle Hamiltonian and the results are 
contained in the fifth chapter. The final chapter concerns a harmonic po- 
tential Hamiltonian. In both cases, the covariance matrices are determined 
and the logarithmic negativity obtained. This provides some insight into the 
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influence of a harmonic potential over the sudden death of entanglement 
Finally, this work's main results are summarised in a short conclusion. 
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Gaussian States 



The Von Neumann entropy is an essential tool in the matter of study- 
ing the entanglement of pure states. In the current chapter, the entropy 
is proved invariant under closed system dynamics, showing that pure state 
entanglement is conserved in closed systems. A formalism is also introduced 
that allows for an easy study of Gaussian states. 



2.1 The entanglement entropy of a two-particle 
Gaussian state 

Let us consider the Gaussian state for two particles in one dimension, 



suggested by Ford and O'Connell FO08 and given by 



_ , . 1 (3=l-^2) 2 (zi+z?) 2 . . 

jff(x 1 ,x 2 ) = e S 2 (2.1) 

v 2Tcsd 



The corresponding density matrix is 

X 2 ] X[, X' 2 ) = H e -^ 2 +-2 2 +x' 1 2 +x' 2 2 )+2e.(x lX2+X ' lX ' 2 ) ^ ^ 
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where Q = e+ = ^ + 1^2 and e_ = ^ - ^ The reduced density 

matrix is obtained by tracing over the position of particle 2 

p 1 (x 1 ,x' 1 ) = J p(x 1 ,x 2 ;x' 1 ,x 2 )dx2 (2.3) 

Explicit calculations give 

p 1 (x l ,x' 1 ) = tte-^ x * +x '^ J e- 2e+ ^ +2e - (xi+x '^ dx 2 

e 2 e 2 I 

with w = e + — z/ = — and fi' = fl./^-. 



2e + ' e + Y 2e 4 

To calculate the Von Neumann entropy, one must calculate the eigenval- 
ues for the state. The general eigenvalue equation is 

J p 1 (x,y)${y)dy = \$(x) (2.5) 

Let us try ^(x[) = e - ^' 1 as an eigenvector. 

J p 1 (x,x')e-< x ' 2 dx' = Q'e^ x2 J e -(«"K)* ,a +^' ' dx' 



IT __„2 , v 2 x 2 



fl\ — — e~ mx (2.6) 
' vj + q 



For this $ to be an eigenfunction, we must have 



v 2 



4(G7 + q) 
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and therefore 



Hence we have the eigenvalue 



Ao = fi\/^— (2.9) 



It is likely that the other eigenfunctions are given by Hermite polynomials, 
which can be defined by 



H n (x) = (-l)V* ( -ft) e~ x2 (2.10) 



For example, H (x) = 1, H\{x) = 2x, H 2 (x) = 4x 2 — 2, H 3 (x) = 8s 3 — 12a;, 
etc. The Hermite polynomials correspond to the Wick objects defined in 



Bar 74 for the arbitrary variable / as follows. 



:/°: = 1 

§- f ■ r . = » . r • : 

<:f n :> = (2.11) 
for n = 0,1,2, .... As a simple example, let us consider two successive terms. 



and 



jj:f:=2(f-(f)) 
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so by integration 



:f := f 2 -2(f)f-(f)+2(fY 



We have the following relation : 



:f n :=c n H n 



f_ 

2c 



where c 2 =< f 2 >. The generating function for these objects is : 



zf 



n=0 



n 



! < e*f > ' 



where 



< e 1 >= exp 
Writing $(y) = e~ ?y2 , we have 



2 ' 2 < / 2 > 



n=0 



Z 



: y n : $(y)dy 

=n'J2~ -y n ■ e -™ x2 -™ y2+uxy - qy2 d y 

n=0 

=tt'e-^ z2<y2> e-™ x2 J e zy e-^ y2+vxy ^ y2 dy 



=tfe-^ z<y> e 



=n'e—* z<v> e-^ J - 



w + s 



-g 4(n7+c) 



-g4(ro+i) g 2(ro+s) g 4(ro 



--Q'J n e ~^ z2<y2> e*(™+^e 2 (™+^$(x), 
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From the Wick objects' generating function : 



cxp 



vxz 



2{w + q) 



-: cxp 



vxz 



:< exp 



vxz 



2(w + q) ' ^ 2(w + q) 



> 



(2.16) 



with 



< exp 



vxz 



2(w + q) 

Hence the previous result becomes : 



>= exp 



z 2 v 2 < x 2 > 
8(w + q) 2 



(2.17) 



oo „ 

/ Pi(x,y) : y 11 : $(y)dy 

n=0 J 

--ft' J— — e~^ z2<x2> 'e _sxi ''e^+^ : exp 



vxz 



zu + q 



2(w + q) ' 



: exp 



z 2 v 2 < x 2 > 
8(w + 2 



7T 



9 2 2^2^ 



12 2 z X — 

$(a;)e~2 z <x > e^+«)e st^+o 2 \ 



n=0 



2(tu + <;)y n! 



: : 



(2.18) 



This becomes of the form Y^=o^n^ '■ x n : W(x), i.e. the z-dependence 
outside the sum cancels provided that 



i .e. if 



< * 2 >= n(m +v-w < 2 - 20 > 
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The eigenvalues are then 



Ar, — £Ia 



7T 2 / U/2 

2e + (w + q) \w + q 



(2.21) 



Using w + q = \Jw 2 — u 2 /4 + w = + v/2 + s/w — v/2) 2 and w + q 



v/2 = y?w -v/2(^/w + v/2 + y/w-v/2), we can check that J2 n = 1 



as 



n=0 



7T 



ra=0 



z//2 V 



7T 



1 



zu + q i 



v/2 



7T 



7T 



+ <T 



2e + V + ^ + — v/2 



vu + q 



zu + q-v/2 2 



7T 



2tt 



zu + q , — . / ? 

+ V e + 



ft'- 1 



2 V~ ' e\- el 



2n 



e+ + \ e+ 



= 1 



(2.22) 



(2.23) 



To calculate the entropy, let us introduce a lemma. 

Lemma 1 Suppose that the eigenvalues of a density matrix p are given by 
A n = £p n where e = 1 — p by normalisation. Then the Von Neumann entropy 
is given by 

oo 

S(p) = - y A n In X n = - ln(l - p) - In p (2.24) 
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Proof: It is quite easy to see that 



S(p) = — A n In X n = — In e — e In p np n 

lne — elnp- 



n=0 n=0 

P 



;i-p) 2 

-lne-e— — hip (2.25) 

1 -pi - p 



One can recall that 



oo 1 

?1=0 



to notice 



5> = T^ = 1 
/ — i 1 — p 

Hence one can write e = 1 — p and conclude the proof. □ 
The entanglement entropy for the initial state is then given by 

S{ Pl ) = - f; A n In A n = - In ( 1 - -^—] - -SEL- In f 

(2.26) 



The entropy is plotted in Figure 2.1 



As expected from the original state, the entanglement entropy vanishes 
for s = 2d and increases again as the separation s increases. 

Proposition 1 The Von Neumann entropy is invariant under closed system 
dynamics. 
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Proof: As a short proof, let us recall that since 



pit) = e iHt Po e-* m 

we have 

S(t) = Tr [-p(t) ln(p(f))] 

Let us call f(p(t)) = p(t)ln(p(t)) and apply the Weierstrass approximation 
theorem |Kre06| so that we can write, in terms of an arbitrary smallness 
parameter e, 



k 



k=0 

p(tf = e tm p e- im e lHt p e- lHt ...e lHt p e- lHt 



which as e — > becomes f(p(t)) = e f(p )e . Then we can write 

S(t) =Tr[-p(t)ln(p(t))] 

= -Tr[e lHt p Hp )e- lHt } 
S(t)=-Tr[p ln(p Q )} (2.27) 

which concludes this proof. □ 

The Von Neumann entropy is well suited to pure entangled states but 
does not suffice to quantify the entanglement of mixed states, such as those 



obtained through open system dynamics RK97 . The logarithmic negativity 
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is a common choice for Gaussian states amongst the many entanglement 



measures that have been proposed KB05,Woo98,MW02 because it is easy to 
compute, especially in the covariance matrix formalism that will be presented 
in the next section. 



2.2 General Gaussian states 

Continuous variables are of growing interest in the field of quantum 
optics and Gaussian states are probably the most widely used states for such 
variables. Indeed, while theoretically easy to handle, the latter can also be 
experimentally prepared and manipulated. A common experimental repre- 
sentation of Gaussian states would be modes of light. Another major advan- 
tage of Gaussian states is their relatively simple mathematical formulation 
and the fact that this formulation allows for explicit calculations. 

Continuous variables 

A quantum system of N particles, each with one degree of freedom, has 
position and momentum variables Rj satisfying (with h — 1) 



Rj, R k 



ia jk l j,k = 1...2N (2.28) 



where R = (xi,pi, xn,Pn) t , x and p are the usual canonical position and 
momentum and a is the fixed, non-singular, skew symmetric matrix defined 
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as 



N 
J'=l 



' o ^ 



One can note two additional properties of a, namely det(cr) = 1 and a 1 



T 

a 1 = —a. 



The state of the system is best described by its density matrix. The 
first moments are collected in the displacement vector d as 



d . = Tr[pi?,] = (Rj) 



(2.29) 



The second moments are collected in the 2N x 2N real, symmetric covariance 
matrix 7 through 



Ijk = Tr 



(Rj - (Rj)), (Rk - (Rk)) 



Ijk + 10 jk = 2Tr[p(Rj - (Rj))(R k - (R k ))} 



(2.30) 
(2.31) 



(Note that the imaginary part is equal to a by the commutation relations 
( |2.28[ )). Not all such matrices are proper covariance matrices ; by (2.30), 
they must in addition satisfy r ) + i<j > 0, i.e. 7 + 20" must be positive definite. 
The covariance matrix 7 is in particular given by 



7jfc = 2Re Tr pfe - {Rj)){R k - (R k )) 



(2.32) 



It may be useful to introduce the Weyl operators W^^ = e l ^ TaR , in terms 
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of which a state p of N modes can be expressed as 

Pr = 7^w[ d2N ^x(-0W^ (2.33) 

where £ is a vector over the phase space R 2N . x( — £) is called the character- 
istic function and is defined as x(C) = Tr pW^ 

Symplectic transformations 

A transformation S on a quantum mechanical state is called symplectic 
if it leaves the canonical commutation relations unchanged. So if 

S : R B! = SR 

is a real linear transformation such that [i?'-, R' k ~\ = tajkt, then SaS T = a. 

Theorem 2 Any real, symmetric, positive definite matrix A can be trans- 
formed into its diagonal form (the so-called Williamson normal form) via a 
symplectic transformation S 

A WF = SAS T = diag(a 1; a\, a N , a N ) (2.34) 

where the a- 's are the symplectic eigenvalues of A. 



The proof of this theorem can be found in Wil36| . The symplectic 



eigenvalues a ? - can be calculated as the positive eigenvalues of ioA. In fact, 
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using a 1 = —a 



a 



1 = S'^a-'S' 1 



S T a = aS- 1 



(2.35) 



eigiiaA) =eig(zaS l A WF S lT ) = eigiiS 1 aA WF S 1-1 ) = eig{iaA WF ) 



-i T \ 



(2.36) 



We can also write 



AV X =XV X 



A 2 V\ =AAV X = \AV X = X 2 V X 



(2.37) 



enabling us to see that the a/s can also be easily calculated as the positive 
square root of the matrix —crAcrA. 



Gaussian States 

A Gaussian state p with N modes is a state whose characteristic func- 
tion Xp(0 can be written as 



X P (0 = exp 



1 

t 



(2.38) 
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where T = a^a T is the covariance matrix of the state and D = ad the 
displacement vector. Considering 



S(j + ia)S T = S 1 S T + icr > 



(2.39) 



one can write the Heisenberg uncertainty relation for r as T + 10 > 0. It 
follows that it can be brought to Williamson normal form by symplectic 
transformation without the symplectic eigenvalues 



Since the matrix must be positive, all the eigenvalues /ifc, k = 1, 2N must 
be positive. They are determined by the usual characteristic polynomial 

N 

J] [(r,- - fi k ) 2 - 1] = (2.40) 

J'=l 

which one can use to determine that 



S[F + ia]S T = Y WF + la = 



i 



j 



(Fj - /i fc ) 2 - 1 =0 



Tj - Hk = ± 1 



/i fc =Tj ± 1 > VA; = 1...2N 



(2.41) 



Hence the symplectic eigenvalues of T are 



r, >i vj = i,...,N 



(2.42) 
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Thus the symplectic eigenvalues of a T satisfying the Heisenberg uncer- 
tainty relations are greater or equal to one. Since the symplectic eigenvalues 
of T can be calculated as the positive square root of Aj where Aj are the 
eigenvalues of —<j r ja r y ) then the Aj > 1. 

Another important property is that since the first moments of a Gaussian 
state can always be made to vanish via local operations, they are irrelevant 
in the determination of entanglement. Hence for the remainder of the present 
work, covariance matrices will be expressed as follows 



7 ifc = 2Re Tr 



pRjRk 



(2.43) 



It should also be noted that due to conflicting notation, the letter 7 will in 
the subsequent chapters refer solely to the coupling constant, whereas G will 
be used to designate the covariance matrix. 

Separability 



We recall (1.8) to determine if a state is separable. In terms of covari- 



ance matrices, separability of Gaussian states can be written as 



Theorem 3 A Gaussian state with covariance matrix 7 is separable if there 
exists covariance matrices 71 and 72 such that 



7> 



71 



V 



72 



(2.44) 
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A proof of this theorem can be found in And03 . For low-dimensional 
systems, one can also check the positivity of the partial transpose (PPT). 
Partial transposition effectively results in time reversal for one of the parti- 
cles'momentum operators. Hence if we have a bipartite system with operators 
Xi,Pi,X 2 ,p2, partial transposition over the first particle results in sending 
Pi to —Pi. With this, one can write 



And03 



Theorem 4 Let T Ta be the partially transposed covariance matrix of a state 
p. I/T Ta fails to fulfill the Heisenberg uncertainty relation, then the state is 
entangled. 

As was seen in the previous section, satisfying the Heisenberg uncer- 
tainty relation is equivalent to restricting the symplectic eigenvalues to values 
greater than one. Hence the above theorem can written a follows. 

Theorem 5 Let T Ta be the partially transposed covariance matrix of a state 
p. If one of its symplectic eigenvalues is less than 1, then the state is entan- 
gled. 



For more formal definitions, one should refer to And03, WW01, SimOO 



VW02]. Although the PPT criterion is not sufficient to completely establish 



entanglement when the system is of dimensions greater than 3x3 WW01 
these conditions allow us to determine qualitatively whether the state studied 
is entangled. To quantify the degree of entanglement in a system, one requires 
an entanglement measure. 

If we write a Gaussian state in terms of its covariance matrix and limit 
our study to second moments, we have 7^ = 2Re Ti[pRjRk\. Partial trans- 
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position may result in one or more eigenvalues to be less than 1, so that the 
resulting covariance matrix may not be positive. Yet, it can still be brought 



to Williamson normal form VW02 , so that its symplectic eigenvalues Xf can 
be calculated using — <j r y Tl <j r j Tl . The logarithmic negativity is then defined 
by 



n+m 



^(P) = -X>g 2 (min(l,| \? I)) (2.45) 

i=l 

This expression makes Ctf(p) easy to calculate and will be used for the 
remainder of the present work. It is again easy to see that since for a separable 
state Xj 1 > 1, Cjsfip) = 0. 
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Entanglement entropy at t=0 
5 | . . . . . . 



4.5 



4 




s 



Figure 2.1: Entanglement entropy for the state at t — and with d = 2 



The Master Equation 



In this chapter, the Non-Rotating Wave master equation is derived 
for a single particle system. The method itself follows that of Gardiner et 



al. MG96GZ00 , but uses perturbation theory as in L096 as opposed to 



the Van Kampen cumulant expansion Kam82 . It starts by deriving the 



Quantum Langevin Equation as done in |LQ88, MG96, GZ00| . The adjoint 
equation follows and finally the perturbation method is used to arrive at the 
master equation expressed in the Non-Rotating Wave approximation and 
quantum Brownian limit. The second part of this chapter gives a solution to 
the master equation for a free-particle system Hamiltonian. 



3.1 Derivation of the Master Equation 



The Quantum Langevin Equation 



To derive the Quantum Langevin Equation, we follow |L088 . Let us 



consider an independent oscillator heat bath, which we couple to a particle. 
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The resulting Hamiltonian has the following form 

3 

where the sum is over all of the bath's degrees of freedom. The bath operators^- 
and qj satisfy the commutations relations [pj,^] = fe,^] = 0. The Heisen- 

2 

berg equations of motion are (with H s = ^ + V(x)) 



■ 1 \u l P 
x =t [H s ,x\ = — 

n m 

1 [H s ,p] = -V'(x) + Y, m i^M ~ x ) ( 3 - 2 ) 



V 



h 



3 



Pj =j l [H s ,Pj] = -mjU)?(qj - x) (3.3) 



It follows that 

rrijqj = - mjUJ^qj - x) (3.4) 
This equation can be solved for qj(t) in terms of x(t) 

qj {t) =q 1 ]{t) + x(t) - / cos [uj(t - t')} x(t') dt' (3.5) 

J — oo 

q*(t) = qj cos(w j t) + sin(wj*) (3.6) 
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Introducing 



fj,(t) — rrijU? cos (cQjt) Q (t) (3.7) 
j 



and inserting it in (3.2) yields the Quantum Langevin Equation 



t 



mx + //(* - ^' + V'{x) = £(t) (3.9) 

J — oo 

where the dot denotes the derivative with respect to time and the prime that 
with respect to x. fi(t) and £(t) describe the interaction with the system as, 
respectively, a memory function and an operator- valued random force. 0(i) 
is the Heaviside function. 

At t — > — oo, one can assume the bath to be in thermal equilibrium at 
temperature T. Then one can write 



( ^ fe)= 2^- C ° th (^)^ 

Ow*>=-pcou.^j*. 



(QjPk) = ~ (PjQk) = -jSjk (3.10) 



The full derivation of the above relations is recalled in Appendix A.l The au 



tocorrelation of the random force, i.e. the expectation of the anti-commutator 
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can then be written as (the [, ]+ denote anticommutation) 



^(EE^[#)^(01 + > 

3 j' 

-(S~^ m 2 u* ( ( a, cos( uud) H — ^— smiud)) ( q-j/ cos(ujd') H — 



sm(Wj. 



+ ( qy cos(wj/t') H — — sin(ujd') ) ( qj cos(ujt) H ^— sin(wjt) ) V 

V ' / V u j m j ) ' 

= - m | w j ( cos(u;ji) cos^t') + 2 — s i n ( w j^) sin(w : ,t / ) 

3 J J 

H (cos(o;jt) sin^t') + cos^-t') sin^t)) 



H (cos^t') sin(cjjt) + cos{uijt) sm(u}jt')) 



rrijUJj 



-\ E ( 2 ^- (H) «» M* - 01 

+ - — sin + *')] - - — sin [uAt + i')]) 



to get 

^«U(01 + > = ^E^l coth (^|) «*M*-0] (3-ii) 

If one introduces jl, the Fourier trasnform of the memory function as 

POD 

fl(z) = / dte tzt fx(z) 
Jo 
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and the spectral distribution G(u) — 3ft [fx (u + iO + )] as 



one can write the autocorrelation of £(£) as 



£(0L> =" rGM^coth ) oosKt-O] (3.13) 



7T 



2kT 



For an arbitrary observable Y of the small system, the Heisenberg equations 
of motion read 



Y 



% 

h 

h 

% 
h 

h 



3 



[x, Y] , / dt' cos [^(t - if)] x(t) 



i[/r.,y]-itt«,y],{(«)] + + i 



/ fi(t - t')x(t') dt' 



(3.14) 



If we work with a Ohmic heat bath LQ88| , then 



/ n{t')x{t') dt' 7±(t) and G(w) ->• 7 (3.15) 
J —00 
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so that the Quantum Langevin Equation for the observable Y is 



Y =- [H s , Y] - JL [[ x , Y ] , at)} + + F] ' 1±{t)] + (3 ' 16) 



The adjoint equation 

The Quantum Langevin Equation is an equation for the system opera- 
tors, whereas a master equation is an (approximate) equation acting on the 
density operator of the quantum system under study. The adjoint equation 
provides a link between the two formalisms, being an exact equation upon 
which approximations can be made to obtain the required master equation. 
One can write p(t) = v{t)pB where u(t) is the density matrix of the small 
system and pB that of the bath. One defines 

Tr s {Y(t)p} = Tr s {Yp(t)} (3.17) 

where Tr s is the trace over the (small) system and Y(t) is a random system 



observable. One then applies (3.17) to (3.16). Term by term analysis yields 



^Tr s {[H s ,Y(t))p} = l -Ti s {[H s ,Y]p(t)} = -^Tr s {Y [H s , p(t)]} 



(3.18) 



^.Tr s {[[x(t),Y(t)],m} + p} 
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= ^_Tr s {[[x,Y],at)] + P(t)} 

= ^Tr s {xY£(t)p(t) - Yx£(t)p(t) + £(t)xYp(t) - £(t)Yxp(t)} 
= ^s{Y[iat),P(t)] + ,x]} (3.19) 



^Tr s {[[x(t),Y(t)] n x(t)] + P} 
= ^_Tr s {[[x,Y}, 1 x} +P (t)} 

= ^Tr s {xY^xp{t) — Yx^xp{t) + jxxYp(t) — ^xYxp(t)} 
= gTr s {Y[[x, P (t)] + ,x]} (3.20) 



Nota Bene The transition from the first line to the second line in (3.19) 
may require some explaining. The trace is over the system's variables ; how- 
ever, £(£) is a bath operator, since it effectively represents the noise. Hence 
it remains unaffected by the trace and retains its time dependency. 

The adjoint equation can finally be written as 
p(t) = ~ [H.,p(t)] - ^ [[t{t),p{t)] + ,x] + ^ [[x,p(t)] + ,x] (3.21) 



The master equation 

To derive the master equation from the adjoint equation, we follow 
[LQ96] . The noise is assumed to be small. This assumption is not essential 
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but merely allows for a simple derivation. To indicate this, we temporarily 
introduce a small parameter e and replace — > e£(t). One can then write 
p(t) to second order in e as 

p(t)=p (t) + m (t) + e 2 p 2 (t) 



The bath and the system are also assumed to be initially decoupled at 
t — > — oo and the bath to be large so that it remains at thermal equilibrium 
throughout. Thus one can write po(t) = uq {t)pB, where ps is the equilib- 
rium state of the bath at temperature T. This assumption is essential to the 
derivation. Then 

p{t)=p,{t) + tp 1 {t) + t 2 p 2 {t) 



The expansion of (3.21) yields 



P(t)=-- h lHs,Po(t))-^e 



[[£(t),p (t)] + ,x] + ^[[x,p (t)} + ,x] 



t[H 9 ,Pi(t)] 



2h 



[m),Pi(t)] + ,x}+^_e[[x, Pl (t)] + ,x} 
[[x,p 2 (t)] + ,x] (3.22) 



Since p (t) = iso(t)pB, (3.22) becomes (with reordered terms) 



% try 

p{t) = ~ ft [Hs, Mt)\ PB + 7^ [[X, Vo(t)] + ,X\ pB 



M t )PB} + , x] + — e [[x, pi(t)} + , x] 



--e 2 [H s , P2 (t)} 



2h 



[[£(*), Pi(t)] + ,z] + 



«T 2 

2h e 



[[x,p 2 (t)] + ,x] 

(3.23) 
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Let us study each order separately. The zeroth order yields 

po(t) = u (t)p B -> Tr B {p (t)} = (t) (3.24) 

so that v (t) satisfies the following equation 

u (t) = - l - [H s ,u (t)] + g [[x,u (t)] + ,x] (3.25) 

The first order term is 

Pi(t) = -\ [H s , Pl (t)} - ± [[£(t)Mt)PB] + ,x] + g [[x, Pl (t)] + ,x] (3.26) 

Since u is a system operator and £ is a bath operator, one can rewrite 
[[£(*)> *4>(*)Pb]+, as 



[[C(*),^o(*)PB] + ,a;] = [K(*),PB]+fo,a;] = [CW,Pb]+NW, 



j; 



so that the first order term becomes 



Pi(t) = -t [H s ,Pi(t)\ + - - m,PB] + [Mt),x] (3.27) 



This can be written in the form 



ptit) = A sPl (t) + f(t) (3.28) 



where 

i 

2h 



f(t) = --\i(t),p B } + [v (t),x] (3.29) 
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and A s is a "super operator" given by 



A sP i(t) = -- [H 9 , Pl (t)] + ^ [[x :Pl (t)] + ,x] 



2h 



(3.30) 



Since Tr# {£(i)ps} — 0, it can be seen that Tig {pi{t)} satisfies the same 
equation as Vo(t), but with vanishing initial condition. It is then easy to see 
that p\(t) can be written as 



Pi(t)= e^WfWtf 



2h 



,Mt-t>) r /,/ 



[voW),x\W),Pb\+ dt (3.31) 



The second order term is 



hit) = [H s , p 2 (t)) - JL Pl (t)) + ,x]+%: [[x, p 2 {t)} + , x] (3.32) 



Using (3.31) 



[[£(*), Pi(f)] + ,s] =- 



I 

2~h 



— oo 

t 



Z(t),e A °W[v (t'),x}[!i(t'),p B U 
[mMt'),PB] + } + e A °W [Mt%x],x 



dt' 
dt' 
dt' 



(3.33) 
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Inserting into (3.32) yields 



p 2 (t) [H s ,p 2 (t)} + ^ [[x,p 2 (t)] + ,x] 



Ah 2 



[Mt'),x],x] [£(*), im,PB) + ] + dt' (3.34) 



Taking the trace over the bath 



h 



Tr B {[H S7 p 2 (t)]} + gTr B { [[±, p 2 (t)] + , x] } 



Ah 2 



Tr 



B 



Ah 2 



[u (t%x},x\Tr B {[C{t),[^t%PB) + ]_ 



dt' 

(3.35) 



Since 



mMt% pb] + ] + \ =2Tr B {m,m] + pb} =2(m,m] + ) 

(3.36) 



(3.35) becomes 



v 2 (t) =-- [H s , u 2 {t)] + ^ [[x, u 2 (t)\ + , x] 

- 1 f [e A ^ [Mt'),x], x ] lam.mu) (3.37) 
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Finally, it remains to recombine all the terms together 



u(t) =Tr B {p(t)} 

u(t) =Tr B {p (t)} + eTi B {^(t)} + e 2 Tr B {p 2 (t)} 
=z> (t) +evi(t) + e 2 z> 2 (t) 

= - 1 - [H s , u (t)] - % - e [H a , v x {t)) - l -e 2 [H s , u 2 {t)] 



+ ^ [[x,Mt)}+,x] + [[x,vi (<)]+, x] 



+ ^e 2 [[x,P 2 (t)] + ,x] 



e 

V 2 



-j_[H s ,v(t)] + ^_[[x,v(t)] + ,x] 



If we recall our assumption that the bath be Ohmic and have a closer look 



at the autocorrelation function, we can see L088 



cos [u(t — if)] dco 



7 



7T 



fto; coth 



2A?T 



cos [u(t — t')] duj 



, ^ d ( TrkT , .. 

=*^ s «*i.(-jj-(*-0 



(3.38) 



which in the classical limit (h — > 0) IFO96 FO02 becomes 



-M(tU(t')] + )^2kTj6(t-t') 



(3.39) 
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Finally, inserting (3.39) into (3.38) yields the Non- Rotating Wave master 
equation in the Brownian motion limit (removing e) 



u{t) = -1 [H s , u(t)) + g [[x, u{t)] + ,x]- k ^ L Mt), x) , x) (3.40) 



For a two-particle system, with each particle in its own heat bath, (3.40) 
reads analogously 

Kt) = ~\ [H s , v{t)\ + ^ [[xt, u(t)} + , Xl ] + ^ [[x 2 , v(t)} + , x 2 ] 

[[v(t),x 1 ],x 1 )- k ^[[v(t),x 2 ),x 2 ] (3.41) 



ft 2 LL'VV.-xj,^ ^ 

where we assume that though the baths follow the same dynamics, they are 
independent and are not necessarily at the same temperature nor have the 
same coupling constant 7. 

3.2 Solution to the Master Equation 

Let us recall the N.R.W. master equation for a single particle in the 



Brownian motion limit {y — >■ p), rearranging (3.40) 



p(t) = -\ [H sys , p(t)) - l g [x, [x, p(t)) + ] - ^ [x, [x, p(t))) (3.42) 



We consider the case where the particle undergoes free motion H sys = 2— 
(and the coupling is a position coupling X x). Considering the position 
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matrix (x\p\y) and p\x) = —ih-^x) (3.42) reads 



d ih ( d 2 d 2 \ 

dt {xlply)= 2^W 2 -w) {xlply) 

_ J_, x _ \(_^__^_ 
2m \ dx dy 



{x\p\y} - -&-(x - y) 2 (x\p\y) 



h 2 



(3.43) 



Applying the change of variables x = u + hz, y = u — hz, (x\p\y) — > P(u, z) 
yields 



d _ 1 / d Id 
dx 2 \ du hdz 



d _ 1 / d _ 1 d 
dy 2 \du hdz 



and 



d_ 

dl 



P{u,z,t) 



d 



2m [dudz^j m dz ^ 



P{u,z,t) (3.44) 



We now apply a Fourier transform with respect to u 



P(g, z ) = / due- iqu P(u, z) 



(3.45) 



to get 



— 2 + — J — + 47fcT2; 
m 2m J dz 



P{q,z,t) (3.46) 



This equation can be solved using the method of characteristics Far82 . The 
characteristic equation is 

dz = 7 z q 
dt m 2m 
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Writing 



we have 



which can be solved as 



w = z + 



27 



dw 7 
— = —w 
at m 



w = ce 



ft/m 



On a characteristic 



d ~ 



dt 



P(q,z(t),t) =-AkT'yz(t)P 



-4ryhT [ce^ t/m - 



P(q,z,t) 

(3.47) 



where c = z + ^. Then 



/ 



dP 



= - 4>ykT 



7 47^ 



In P - In P = - AjkT 
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'y 2 47 2 



P(q,z(t),t) =P(q,z ,0)exp 
x exp 



/cTt 



7 



■5 



-2mkTc 2 (e 2 ^ m - 1) + eg (e^ m - 1) 

7 
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To obtain P(q, z, t), we now have to express c in terms of z and t with 



<■ = I z - * ) e-' 1 



Finally the solution can be written as 



P(q,z,t) =P(q ,z ,0)exp 



with 



kTt 



+ 



7 

AmkT 



q 2 - 2mkT (z + (1 - e~^ t/m ) 



7 



: • £ ) q{\ , 



--P(q Q , z Q , 0) exp 



- rtq 1 



(3.48) 



(3 =2mkT(l - e - 2lt/m ) 

_Q_ 
27 



7 

zo =ze~ 

kT 

7 



r 



(3.49) 
(3.50) 

(3.51) 

(3.52) 



We have thus found a simple solution to the non-rotating wave master 
equation, in the case that the system Hamiltonian is chosen to be a free 
particle one. Then we need only apply a Fourier transform to the initial 



state and substitute for P(qo, Zq,0) in (3.48). The method used to derive 



(3.48) can also be used in the case where a harmonic potential is added to 
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the system Hamiltonian, as will be shown in Chapter 4 in the one-particle 
setting and in the final chapter for a two-particle case. 

3.3 Comments on the coupling constant 

This short section is intended to give a few lines of comments on the 
coupling constant 7. As shall be seen in the following chapter, the bigger 7 
is, the stronger the coupling. One may, for instance in the case of a harmonic 
oscillator with frequency oj , notice two different damping behaviours (under 
and over-damped), depending on how large 7 is with respect to ojq. The 
reservoir considered is always assumed to be large compared to the system 
and as such is not affected by the system's influence, i.e. one can assume 
that the bath remains essentially in thermal equilibrium (since we are using 
thermal baths) and as such its state is roughly time-independent. This means 
that at any time we may write p(t) ~ p s (t) <g> ps- However, the evolution 
of the system will in general be strongly influenced by the coupling to the 
reservoir. Thus strong coupling is to be understood as strongly affecting the 
system. The noise is still small in the sense that it affects the reservoir little. 
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The Non-Rotating Wave Master 
Equation for a Single Particle 



This chapter studies the time evolution of a single-particle Gaussian 
state, comparing the results obtained with the N.R.W. master equation to 
the results obtained by Savage and Walls in [SW85b| . It also examines the 



case of a quantum harmonic oscillator and compares to work done by Savage 



and Walls in SW85a|. This will allow us to determine that the master 



equation chosen is satisfactory to our purpose. 

4.1 The case of a free particle coupled to a heat 
bath 

Consider a one-particle system with a Gaussian state wavefunction. 
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The corresponding density matrix is 



'7TS Z 



(4.2) 



We evolve it according to (3.42), with a free particle Hamiltonian. Re- 



name x = u + hz and y = u — Kz so that 



p(x,y) P(u,z) 



TCS Z 



(4.3) 



Applying the Fourier transform 



P(q,z) = / due- iqu P{u,z), 



(4.4) 



yields 



P(q ,z Q ) = exp 



h 2 z 2 qls 2 



(4.5) 



The solution to the master equation is recalled here (see (3.48) 



P(q,z,t) = P(q, z ,0)exp 



rtq 1 



(4.6) 



with 



f3 =2mkT(l - e - 2lt/m ) 



7 

Zq = ze -^ t/m - — (1 - e" 7 * /m ) 
27 v ' 



T 



kT 

7 



(4.7) 
(4.8) 

(4.9) 

(4.10) 
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Substituting into (3.48), we get 



P{q, z, t) = e -^-m- 2 +~s(t)z q ( 4 U n 



with 



s 2 kTt mkT . _„,v-/™n ^ 2 



e 



-7*/m\2 



a{t) =- + —HI - e-^ m )(3 - e'^ m ) + - 

w 4 7 27 2 v >y ' 47 2 s 2 

(4.12) 

=2mkT{\ - e -^ t/m ) + — e" 27i/m (4.13) 



s 



h 2 9mIcT 

S(t) =— e -^ m (l - e~^ m ) + _ e -7*A»)2 (4.i4) 



7s 2 7 



Performing the inverse Fourier transform, we get 



p(x,y,t) 



x exp 



1 nr 



27ry 5(*) 

5(t) 2 \(x-y) 2 {x + yf ,>(/) 



0(f) 



45(f) 



16a(t) 8d(t) 



(4.15) 



Savage and Walls considered the following master equation in the zero 



frequency limit and at high temperature SW85b,SW85a 



d t p = —ih 1 



2m' 



[r, Pp + pP] - 2m-fk B Th- 2 [r, [r, p]] (4.16) 



Using the characteristic equation approach, they solve the equation for 
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the initial state 



tf(r) = (7ra/2)- 3/4 e- r2/ff 



(4.17) 



to get 



(r|p|r -n) = (2nc)~ 3/2 exp [-ibn . (r - /x/2)/c - (r - ^/2) 2 /2c - (a - 6 2 /2c)^ 2 ] 

(4.18) 

with 



mkT 



2h 2 
kT 



1 - e" 47 ' 



ft 2 



mkTa 



2-fh 
a 

4 + 



(1 



,-27* 



) 1 



-27* 



/i 2 



ft 

2m7 



+ 



mkTa 
kT 



a 



mj^ 



7t- 



,-27* 



-47* 



(4.19) 



To compare (4.12) and (4.19), one may first remark the difference in 



the definition of the momentum term of (3.42) and (4.16). This leads to 
e -7t/m ^ e~ 2lt . Then if one notes that a — » s 2 , one can see easily that c ~ a, 
a ~ (3 and 6 ~ 5. This allows us to say that both master equations are in 
exact correspondence. 
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4.2 Study of a one-particle Gaussian state with 
a harmonic potential 



Although as we have seen, the master equation seems to be consistent 
with past results, it is useful to verify it for a particular case. The dynamics 



of the harmonic oscillator are very well known Sen60 Sen61 ,HW85 and thus 
make it a perfect candidate for study. 



Evolving the density matrix 



Let us evolve (4.2) according to (3.42), where the harmonic potential is de- 



fined as H sys = y~ + oox 2 . The master equation then becomes 



?E---\rr 1 _ 
dt~ h [sys,Pl 2mh 



[x, [p,p] + ] - -p- [x, [x,p]} 



(4.20) 



Using the position space matrix form (x\p\y) of the density matrix, (4.20) 
reads 



0_ 
Of 



\Ap\v) 

ih d 2 d 2 ioj 
^dx 2 ~ dy^' ~ 2h 

7 d d 

— (x — Z/) ( — — — 
2m dx dy ' 



f I A 



x-yf \{x\p\y) 



(4.21) 
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which we now solve. Rename x = u + Hz and y = u — hz so that (4.21) 
becomes: 



at \ 2m ouoz m oz J 



Now we write P in terms of its Fourier transform P with respect to u, 



and replace it into (4.22) to get : 



dP + q + 2jz dP _ 2ujz dP = _ 4ikTz 2p 
dt 2m dz dq 



(4.23) 



The characteristic equations are 



a M 

OtV = — V 

2m 



(4.24) 



with v T = (z, q) and 



M 



\ 



2 7 1 
-Amu 



To solve the differential equation, we need the eigenvalues and eigenvec- 
tors of M. Those can readily be calculated : 



X 1 



7 + a/7 2 — 4wm, 7 — a/7 2 — 4um J = (A + , A_) 



(4.25) 



and 



/ 



i 

" A_ A 
1 



>+ 

1 J 
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Since Q 1 MQ = D, we must calculate Q 1 to be 



/ 



A_ - A+ 



A_A+ A_ 
-A_A+ -A+ 



/ 



The differential equation then becomes 



One can then write 



2md t w =Mv = QDQ- 1 v 

a D 
o t w — — w 

2m 
w =Q- 1 v 



v(t) = Qe Dt / 2m Q-\ 



and also 



v = Qe- Dt / 2m Q-\{t) 



with 



This yields 



±Dt/2m 



( e ±\+t/2m q \ 



/ 



v(f) = 



\ +e X + t / 2m — \_ e X - t / 2m e \ + t/2m_ e \_t/2m 



x ) ^(e A +*/ 2m -e A -'/ 2m )^o + 



\_ e x + t / 2m — \ +e x - t / 2m 

XI^aT 9o 



52 



and similarly 



v 



/ \ + e~ X + t/2m -\-e~ X - t/2m e -\ + t/2m_ e -\_t/2m 

A_-A + Z A_-A + 1 

A-A + , e -\+t/2m _ e -\.t/2m) z + X-e^+^-X+e^-^ g 



A_-A+ 



/ 



(4.32) 



The equation for P becomes 



J? 



-4rykTz 2 



4 7 £;T 



x 



( A+e A +'/ 2m -A_e A - t / 2m )% 2 



dP 
P 



+ ( 



3 A + t/2m _ A_t/2m\ 2 2 



\^ 2 
) % 



+ 2z q (A + e A +*/ 2m - A_e A -*/ 2m ) (e A +'/ 2m - e A -'/ 2m ) \p 



(A- - A+) 2 

x | (X 2 + e x ^ m + A 2 e A -'/ m - 2A + A_ e ( A ++ A -^ 2m ) z 2 

+ ^A+t/m + e A_t/™ _ 2e (A++A-)t/2m) ^2 

+ 2^ g (A+e A +*/ m + A_e A -*/ m - (A+ + A_)e( A + +A -)'/ 2m ) \dt 



Then 



lnP-lnP = 



A^mkT 
(A- - A + )- 



x 



^A+e A +' /m + A_e 



A ///„ _ 4A + A_ c (x ++ \_) t /2m 



+ 9o 1 Tf 63 



1 

+ A7 e - 



A+ + A- 
4 



A + + A_ 



3 (A++A_)t/2m 
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P = Pq exp < 



+ 2z q (e x+t/m + e x ~ t/m - 2e {x++x - )t/2m ) 
A'ymkT 



(A- - K) 2 

z 2 (\ + (e x+t/m - 1) + A_(e A "' /m - 1) 



x 



4A+A_ 
A + + A_ 



^ e (\ + +X.)t/2m _ -QA 



X+t/m _ ^ e X-t/m _ -|_ 



(A++A_)t/2m 



+ 



1) 



A_ 



A + + A_ 



+ 2z q (e x+t/m + e x ~ t/m - 2e {x++x - )t/2m ) 



(4.33) 



Let us write (4.32) in the form 



zq = ai(t)z + h(t)q 



qo = a 2 (t)z + b 2 (t)q 



Inserting into (4.5), we have 



P(q , z ) = exp 



ft s 

2zq{—a 1 {t)b x {t) + S -a 2 {t)b 2 {t)) 
s z 4 



(4.34) 



Thus (4.33) becomes 



P = exp [-A{t)q 2 - B(t)z 2 - C(t)zq] 
with 
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Transforming back to the variables x,y, we get p(x,y,t): 



p(x,y,t) 



1 

2tt 



71 



Ait) 



cxp 



- Bit) 



— i 



Citf \ (x 
AA(t)J 
C(t) x 2 



yf 



y 



Ah 2 

1 (x + y) 2 ' 



1 

2tt 



2A(t) Ah 

tx 2 —ly 2 +2vxy 



4A(t) 



A(t) 



(4.39) 



with 



t= ^ m '^))w 2 + im +l ^W) =i+iri 
B w-S)^-«dw (440) 



and e the complex conjugate of e. To check that the density matrix is 
normalised : 



J p(x, x, t)dx = Q J exp (-(e + e - 2v)x 2 ) dx = 1 (4.41) 



if Q = — 

11 " 27rV^(*)- 
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Study of the entropy 

To calculate the entropy, one must get the eigenvalues of the density matrix. 
The eigenvalue equation is the following : 

J p(x,y,tMy)dy = Xo(t)*(x) (4.42) 

Let us try *&(y) = e~ 6y2 as an eigenvector. 
We get 



We must calculate 5, where 



v 2 



5 = e- (4.44) 



We have (5 + e)(5 — e) = — v 2 so 5 2 — 5(e — e) — |e| 2 = — u 2 , i.e. 5 2 + 2t5i] — rf = 
£ 2 — v 2 calling e = £ + irj. 

S(t) = y/l 2 - v 2 + irj (4.45) 

and 

S(t) + t=^e-^ 2 + i- (4.46) 
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The corresponding eigenvalue is given by 



A (t) = n 



6 + e 



(4.47) 



If one recalls (2.9) and (2.21), one can write the other eigenvalues as 



7T 



u{t) 



(t) + 5(t) \e(t) + 5(t) 



(4.48) 



The sum of the eigenvalues must equal 1. In fact, 



n=0 



n=0 



7T 



e + 5e + 5-z/' 



(4.49) 



and using e + 5 = a/£ 2 — z/ 2 + £ = §(\/£ + ^ + vT~~~^) 2 an d e + 5 — v 
V£ ~ v{\fl + v + Vf - we S et 



n=0 



Using Lemma [TJ the entropy is given by 



(4.50) 



S(t) = - In 1 



v 



In 



(4.51) 



Figures (4.1, 4.2) show the entropy in the under-damped and the over- 



damped case respectively, increasing as the equilibrium is destroyed over 
time. One can observe that the over-damped entropy increases more smoothly 
than its under-damped case counterpart. One can also see that the entropy 
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rises to larger values for the over-damped case than for the under-damped 
case. A rather quick analysis of the behaviour of the entropy reveals that at 
large time, e - x + t / 2m — e -^~t/2m = g so that a\ = a 2 = h\ = b 2 = leading to 
A{t) = B(t) = C(t) = regardless of damping. Then 



v 

1 ~ 2 

e + 5 



and 

v AA(t)B(t) - C(tf - h 2 



S + e y/4A(t)B(t) - C{t) 2 + h 
If we recall M's eigenvalues 



1 



7 + a/7 2 — Aojm, 7 — a/7 2 — Aojm ] (4-52) 



we can readily observe the differences in damping. Indeed the bigger 7 is 
with respect to u, the more strongly damped the system is. The eigenvalues 



A + and A_ will be real only in the over-damped case 7 > 4wm. Figures (4.1 



4.2, 4.3, 4.4) illustrate the differences between the types of damping quite 



clearly. One can recall the results of Savage and Walls in SW85a| and study 
the off-diagonal terms by replacing y = x — and getting the variance of the 

2 

e~ M coefficients. Some algebra yields 



Figure 4.3 shows the time evolution of A/i for the under-damped case. 



One can easily notice the decreasing oscillations whereas on Figure 4.4| rep- 
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resenting the over-damped case, one can see that the off-diagonal elements 
disappear quite quickly. Since the off-diagonal elements of the density matrix 
represent coherence |Zur03| , one can easily conclude that the environment 
tends to diagonalise the density matrix, all the faster the stronger the cou- 
pling. This agrees with the result of Savage and Walls S W85 b , SW85a| . 
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Entropy at o= 0.85 and gm = 0.05 




2 4 6 8 10 12 14 16 18 20 

t 



Figure 4.1: Entropy vs t in the highly under-damped case 
This plot is obtained for 7 = 0.05 and u = 0.85 



Entropy at 0= 0.05 and gm = 0.85 

2.5 1 1 1 1 1 1 1 




Off diagonal coefficient at o= 0.85 and gm = 0.05 




2 4 6 8 10 12 14 16 18 20 

t 



Figure 4.3: A/i 2 vs t in the highly under-damped case 
This plot is obtained for 7 = 0.05 and u = 0.85 



Off diagonal coefficient at 0= 0.05 and gm = 0.85 
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Free Particle Dynamics 



The second chapter has established both the conservation of entangle- 
ment in closed systems dynamics and the formalism used for the main aim of 
this work. The present chapter contains one of the main results of this the- 
sis and concerns the evolution of the entanglement in a two-particle system 
when it is left to evolve freely while subjected to an environment. 

5.1 Free particle Hamiltonian 

The time evolution 

Let us recall the entangled Gaussian initial state that we studied in 
Chapter 2 



\Jl-Ksd 



1 



e 



(x 1 -x 2 ) 2 _(x 1 +x 2 ) 2 
il 2 e 16d 2 



(5.1) 



with corresponding density matrix 



p(xi, X2; x'^x^; 0) = Q e 



(x 1 -x 2 ) 2 _ (x x +x 2 ) 2 _(x' 1 -x[ 
16d 2 e is 7 



x 2 ) 2 ( x \~^~ x 2 ) 2 
72 16d 2 



(5.2) 
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The evolution of the state is modelled via the N.R.W. master equation 



Pit) 



- ^ [H„ p] + ^ [[xi,p]+ , xi] + H [[x 2 , p]+ , x 2 ] 



h 2 



(5.3) 



whose solution is recalled here : 



P(q,z,t) = P(q ,z ,0) x exp 

x exp 



where 



2 7l 

272 



-^2(^2 + t^-) 2 + a 2 g2(^2 + ^-) - r 2 tq 2 2 



Qi 
27i 

g2 
272 



(5.4) 



z = 

ft = 

Oti = 

Ti = 



[Zi + 



Qi 



-Jit/m Qi_ 



2%' 2% 
2mkT l (1 - e- 2 ~<> t/m ) 
AmkT; 



li 

m 

it 



1 - e -~ Ht/m ) 



(5.5) 
(5.6) 
(5.7) 

(5.8) 



In order to include the initial state in the solution, we must first perform 
the following change of variable 



x = u + Hz 



x = u — Hz 



(5.9) 
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to get 



P(u, z; 0) = il exp 



(ui + hz 1 — u 2 — hz 2 ) 2 (ui + hz 1 + u 2 + hz 2 f 
_ 

(ui - hz 1 -u 2 + hz 2 ) 2 



16 d 2 

(ui - hzi + u 2 - hz 2 ) 2 



4 s 2 



16 d 2 



(5.10) 



We can simplify the above expression as follows : 



(ui + hzi — u 2 — hz 2 ) 2 + {u\ — hzi — u 2 + hz 2 ) 2 = 2 (ui — u 2 ) 2 + 2 {Kz\ — hz 2 ) 2 
(ui + hzi + u 2 + hz 2 ) 2 + (u\ - hzi + u 2 - hz 2 ) 2 = 2 (m + u 2 ) 2 + 2 (hz 1 + ^ 2 ) 2 

(5.11) 



to get 



P(u, z; 0) = f2 exp 



- (^ + ^)K 2 + + M 2 2 + h 2 z 2 2 ) 

+ 2(- 



•2 s 2 8d 2 



)(lilli2 + ft -21-22) 



(5.12) 



We now apply a Fourier transform to get 
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P(q o ,z o ;0)= J P(u,z ;0) e'^'^dm du 2 

=$7 exp [—e + h 2 Zol — e + h 2 z l + 2e_h 2 z 01 z 02 \ 

x J exp [-e +Ul 2 - e+V + 2e_WiW 2 ] e^ 1 " 1 ^ 92 " 2 du 2 

(5.13) 

where e + = 7^2 + g^j and e_ = ^2 — g^2- The integral with respect to 
u 2 is performed first, yielding 



P(q , z ; 0) =Q exp [-e + h 2 z l - e + H 2 z l + 2e-h 2 z 01 z 02 ] 

-2e_Mi + 2g 2 ) 2 



exp 



-e+ui - iq\U\ + 



4e 



+ 



du\ 



P(q ,z ;0) =fi 



7T 




7T 



exp [—e + h 2 z l — e + h 2 z l + 2e_ft 2 ^01-202] 



exp 



g 2 2 (?i + ^) 2 



4e_ 



4 ( e + - fe) j 



(5.14) 



If we now replace z = (z + £.) e 7 */ m — ^ and using O 
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we get 
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^(q >zo;0) 



exp 
x exp 
x exp 



+ 



4e + ' 4e + (e 2 + - el) J q2 4(e 2 _ - e 2 _)^ 2(4 - e 2 ) 

2 



:9i 



g2 n e - 72 t/m _ 

272 2 72 



2 e _n 2 ((^ + ^)e-^ 
V 2 7i 



gi 
271 



(.2 + ^) e"-^ - ^ 



(5.15) 



Finally 



P(q,z,t) = exp 

x exp 

x exp 

x exp 

x exp 
x exp 



C+ n 2 

-Q2 — 



-q 1 - 



4(4 -ei)^ 4(4 -e 2 )^ 1 2(e 2 - e 2 ) 
_ e+tf ((, 1 + |L )e -^_^ 2 



?1?2 



271 271 



-e + n 2 ((.2 + ^)e-^-^- 
272 272 



-^2(^2 + 7T^-) 2 + "292(^2 + - r 2^2 2 

2 7 2 2 72 

V 2 7l ; 2 7 J V 2 72 ; 2 72 

(5.16) 



We can write this in the simpler form 
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P(q z t) = Q~Aiqi 2 -A2q2 2 -Bizi 2 -B2Z2 2 - r Dz 1 z 2 -£qiq2 

-C\\z\q\— C22Z292— Ci2Zi92— C21Z291 



x e 



with 



^ =77^T + ^ - ^ + A + Pi (1 - ^ 7lt/m ) 2 
4(e 2 - e 2 ) 2 7 i 4 7i 2 4 7i 2 

^ 2 =77^T + ^ - ?- + A + ^ C 1 - e ~ 72t/m ) 2 
4(e 2 - e 2 ) 2 72 4 72 2 4 72 2 

Si =h 2 e + e' 2 ^ t/m + 0! B 2 = h\ + e ~ 2 ^ t/m + 

C u =&- ai - f^U-nt/m _ e -2mt/m^ 

71 7i 

72 " 72 

I? = - 2ft 2 e_ e~ 7l * /m e" 72t/m 

2(4 -e 2 _) 2 7l72 (1 K1 } 

C12 =— e~ 71 '/ m (1 - e~ 72t/m ) C21 = — e" 72 * /m 

7 2 7i 



Applying the inverse transform yields 



P(u,z,f) = 7^/ P(ci,z,t))e^ + ^dq 1 dq 2 
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4tt 2 



-B 1 z1-t3 2 z%-T>z 1 Z2 



,-A\q{-A2qi,-£-q\q2-C\\z\q\-C22Ziq2-C\iz\qi-Ci\Z2q\ 



x exp [iq\Ui + iq 2 u 2 ] dq± dq 2 



1 n e -Bxz\-Bizl-V .:.-._> 



4tt 2 V ^2 

y exp \-A\q\ - (CuZ! + C 21 z 2 - iu ± ) q ± ] 
(C 22 z 2 + C 12 z x + Sq x - iu 2 ) 2 



x exp 



AA 2 



dqi 



IT 



7T 



^]/A 2 



x exp 



x exp 



(c 22 z 2 + c 12 zi - iu 2 y 

4~A 2 



4(A 



£ 2 - 
4.4 2 ' 



4vr 2 



4tt 2 V 4^M 2 -5 2 

C 22 z 2 "I" Cf 2 zf — u\ — 2i (C 22 z 2 + C\ 2 Z\)u 2 + 2 C 22 C\ 2 z 2 Z\ 



x exp 
x exp 



4 (^i - &: 



4 ^ 2 

^11^1 + ^21^2 + 2CnC2l^l^2 



£ 



~ u i + -r^ 2 ( C 22Z 2 - C X2 z x + m 2 y 



AA 2 

£ 



£{C 22 z 2 + C\ 2 z\ - m 2 ) 
2^ 



(C n zi + C 21 z 2 )(C 22 z 2 + Ci2^i - iu 2 ) 
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Finally 



4vr 2 



x exp ( - loi z\ -u 2 z\- xi u\ ~X2u\- X12 u x u 2 - 5 z x z 2 
- i9 n u x z x - %9 22 u 2 z 2 - i9 12 u x z 2 - i9 2X u 2 z x ) 

(5.19) 



with 



A 2 Ai -S 
Xi =-r- A — ^ X2 = -r-7—. ^ X12 



4A X A 2 -S 2 " AA1A2-S 2 ^ 4A X A 2 -S 2 

u x =& x ~ C X2 X2 - C xi xi - C12C11X12 
uj 2 =B 2 — C 2 \x2 — C 2 \xi — C 22 C 2 \X\2 

5 =T> - 2C22C12X2 - 2C11C21X1 - X12 (C n C 22 + C12C21) 
#11 =2CnXi + C12X12 #22 = 2C22X2 + C21X12 

012 =2C 2 iXi + C22X12 ^21 = 2C 12 X2 + C11X12 (5.20) 



Going back to the original variables yields 



— ?ja:ia!2— T/a^aij— C^ia^ - C' 3123 '! 



(5.21) 
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with 



''Ah 2 



+ ^7- +i 



33 

Ah 



Ah 2 



+ -% 



33 

Ah 



Xj 

Ah 2 A 



X12 , 5 Q\i 2 i 
A Ah 2 Ah Ah 



c = 



Xl2 



Ah 2 



&12 , 021 

l Ah +t Ah 



12 







21 



' - Xl2 6 9 

71 ~ A + Ah 2 \h \h 



>/ _ X_12 _ 0\2 021 

~ A Ah 2 Ah 1 Ah 



(5.22) 



Calculating the covariance matrix terms 

We now proceed to compute all the elements of the covariance matrix. In 
all the following, the integral with repect to x 2 is performed first. Note that 



4XiX2 - xi 2 



AAiA 2 -£ 2 



and Q' 



2lT ^4AiA 2 -£ 2 ' 



^ =2Re Tr[pXiXi] 

= 2Re jfi' J x 2 1 e- XlX '- X2X2 - Xl2XlX2 dx 1 dx 2 ^ 
^ = 2Re I Q', 



2X2 



4xiX2 -X? 2 4xiX2 -X? 2 



AAx 
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" d 


| — ih Tr 


OX i 



4?r 2 / 4fa - i/i)x2 _ Xi2(?? + C) \ 

4 Xi X2 ~ X? 2 V 4 XlX2 - X? 2 4 XlX2 - X? 2 / 



=2Re lihQ! 



n. =2Re Tr[pXiX 2 ] 

13 



=2Re ( ff / we ^H^ 



1 = 2 ^ W „ ^ 2 , " Xl2 2 

13 1 V 4xiX2 - X! 2 4 XlX2 - XI 2 



IE 





" d ' 


j-^Tr 


Xl dx 2 P _ 



^ =2Re 



2(^2 - ^)Xl2 2(77 + C')X2 



4 Xi X2 - X? 2 V 4 xiX2 - X? 2 4 xiX2 - X? 2 , 
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=2Re Tr[pAXi] = 2Re j«tfTr[ J 

V 4 xi X2 - X12 V 4 XlX2 - Xl2 4 XlX2 - Xl2 

=2Re Tr[pAA] = 2Re j-ft 2 Tr[£^p]j 

=2Re <fft 2 ft'2^ J e- XlX i- X2X 2- Xl2XlX2 dx 1 dx 2 

h 2 n'A^[ - z/x) 2 J xj e -xixl-X2x 2 -xi2x 1X2dxi dx2 
h 2 n'4(q[ - vJirf + CO J x lX2 e -xi*i-X24-xr2xix 2dxi dx2 



= -c 



=2Re <!n 2 fi' 



4tt 2 



4xiX2 - X? 2 

f 2,' 

V 1 4xiX2 - X? 2 4 XlX2 - X?2 4 XlX2 - X?2 



fQ ^ 8(^-^) 2 X2 4(^-^)^ + 0x12 2(V + C / ) 2 Xi 
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=2Re Tr[pA^ 2 ] = 2Re jzftTr[J^V 2 ]j 

-tMV(rf + C') J x 2 2 e- XlX *- X2X 2- Xl2XlX2 dx 1 dx 2 
[ V 4 Xi X2 - XI 2 V 4 XlX2 - Xl2 4 XlX2 - XI 2 / 

=2Re Tr[pAA] = 2Re {-^T r [AAp]} 

- ft 2 ft'2(4 - i/ 2 )(t/ + C') y r^e"* 1 ^-* 2 ^-* 12 * 1 ^^ 

- h 2 Q'2(q[ - z/x)(V + C) J x 2 1 e- XlX "~ X2X2 - Xl2XlX2 dx 1 dx 2 

- h 2 Q' (4(^ - ^)(<r 2 - 1*) + (1/ + CW + 0) 

■2i\v >h 2 n' 



4tt 2 



x irf - 



4^XiX2- X12 

4( ^-^)^ + 0x2 _ 4(^-^)^ + 0x1 
4xiX2 - X? 2 4 XiX2 - X? 2 

, xi 2 ((V + CW + Q + Wi-viM-i*) ) 

4xiX2 - X? 2 
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-7 =2Re Tr[pX 2 Xi] = ^ =28 

31 13 



^ =2Re Tr[pX 2 A] = 2Re Tr[X 2 Ap] = 2Re \ -ihTr[x 2 ^-]\ 
32 [ ctei J 

3 =2Re ,/ , ^ , f- 2fe -^ 2 + 2 <" + «*' 



32 ! V 4 Xl X2 ~ Xl2 V 4 XlX2 - Xl2 4 XlX2 ~ Xl2 



-1 =2Re Tr[pX 2 X 2 ] 

C 33 



3 



33 



-2Rc UYJ A ^ 2A 2X1 2 \ = U, 

4XlX2 - Xl2 4 XlX2 - Xl2 



^ ^ =2Re Tr[pX 2 P 2 ] = 2Re Tr[X 2 P 2 p] = 2Re |^Tr[x 2 ^]| 

7 =2Re I ^ ~( ^-^ Xu(r, + C) \ 

V 34 1 V 4xiX2-XI 2 V 4 XiX2-XI 2 4xiX2-XI 2 / 
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-7 =2Re Tr[pP 2 ^i] = 2Re IthTi^p^ 

U 41 |^ 

1 =2Re U< ./ , ^ , f ^% - ^±% ) 

41 I V 4 Xl X2 - Xl2 V 4 XlX2 - Xl2 4 XlX2 - Xl2 J 



-l =2Re TV[pAA] =3 = I 

U 42 U 24 Z 



^ ^ =2Re Tr[pP 2 X 2 ] = 2Re j^Tr^-^-] J 

-2Re IMJ . ^ , f- 4 ^~^ + X U W + C) 



^ =2Re J zMT * / — — =- ( 

^43 1 V 4 XlX2-X?2V 



4xiX2-Xl2\ 4xiX2-XI 2 4xiX2-X! 2 , 
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- n 2 fi'4(4 - z/ 2 ) 2 y x 2 e -xi*?-x a ^-xi2xi*2 da . 1 ^ 

-^^(v+o 2 y x 2 e - xi ^- x2:c 2-»2^2 rfxirfX2 

^ 2 fi'W =- 

^ 2 4xiX2-X? 2 4xiX 2 -X?2 4xiX2-X? 2 ^ 



Alternative derivation of the covariance matrix terms 



The covariance matrix terms can also be calculated using the change 



of variables x = u + hz and x' = u — hz, and (5.17) as follows 



{XtXj) =Tr[pX,X J ] 

oc i Xjp(yObi^ oc 'jj doc ^ doc j 



■ J Ui UjP(u, z = 0, t)dui duj 
d d A . 



(5.23) 



if we notice that 



^-P = -i I i,P< 'dii 
oq 



(5.24) 
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Similarly, we can get 



f d 

(X i P j )=-ih J Xi—p(xi,Xj)dxidxj 



1 f d i f d 

ih I Ui— — P(u,z,t)du,iduj / Ui— — P(u, z = 0,t)duiduj 

2 J (ju j 2 J j 

If 1 f dui 

— -%h \ UiP(u,z,t)duiduj + -ih / — —P(u, z,t) 



% f d 

2 J U i-Q^. P ( U ' Z = Q ' t ) du i du i 



If % f d 

= -th5ij / P(u,z,t)duiduj — - / Ui- — P(u,z,t)duiduj 
2 J 2 J OZj 



1 / d d ~ 

^ihbij J P(u,z,t)du i du j + - y——P(q,z = Q,t) ) \ z =o, q =o 



(5.25) 



f 9 

(PiXj) =ih J Xj—p(Xi,Xj)dxidxj 



1 / d d ~ 



- ^th5 i:j J P(u,z = 0,t)duiduj + ^ ^^-^-P(q,z = 0,t) 



1 2=0,9=0 

(5.26) 



and 



d_d_ 

dec 'i dec j 



p^yObi^tj eCj^diCCi dxj 



4 . 

1 ( d d ~ 



d_ ld_ 

dm h dz; 



4 V dzj dzj 



P(q,z,f)^ | 2 = ,g=o 



(5.27) 
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Finally, 



2Re(X i X j ) 
2Re(X i P j ) 
2Re(P i X j ) 
2Re(P i P j ) 



, d d ~, A , 

d d ~. A . 

d d ~, A . 
^% P(q ' Z)t) J 
^ f d d ~ f A . 

(q ' Z '°J 



(5.28) 
(5.29) 
(5.30) 
(5.31) 



This method allows us to get the same results as the explicit calculations 
of the covariance matrix but in a much easier fashion. 



5.2 The Logarithmic Negativity 

Using these, we can write ^ in matrix form : 



1 



4A -Cn 2£ -C 2 i 

-C n fix -Ci 2 P/2 

2£ -C12 4^ 2 -C 22 

— C 2 i "D/2 — C22 ^2 



(5.32) 



To obtain the partial transpose, we set f>\ — > — pi so that ^ b 



ecomes 



1 



4Ai Cn 


2S 


-C21 


C u Bi 


C\2 


-V/2 


28 C12 




—C22 


-C21 -V/2 


—C22 


B 2 



(5.33) 
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In order to use the logarithmic negativity as defined in previous chapters, let 
us calculate a^j T to get 



1 



Cu 


B 1 


C 12 


-V/2 


-AA X 


-Cu 


-28 


C21 


-C21 


-V/2 


—C22 


B 2 


-2E 


—Cyi 


-4:A 2 


C22 



Then we calculate — a 1 ^ T °"~^ T to S et 



-°1 1 = 

where tU 12 = tU 21 = tU 34 = tU 43 = and 



m 


11 


m 


12 


m 


13 


m 


14 


m 


21 


m 


22 


m 


23 


m 


24 


m 


31 


m 


32 


tu 


33 


m 


34 


m 


41 


m 


42 


m 


43 


m 


44 



(5.34) 



(5.35) 



m 11 = m 22 =44iiBi -vs + c l2 c 2 i - c 2 n 
m 33 = m 44 =4A 2 B 2 - c\ 2 -vs + c 12 c 2 i 
m 13 = m 42 =2SB 1 - 2A 2 v - c n c 12 + C12C22 
m 14 = -m 32 = - c 12 B 2 - C21B1 + CuV/2 + C22V/2 

m 23 = -HI 41 = - 2SCu + 4AlC 12 + 4A2C21 - 2£C 2 2 

m 24 = m 31 =2SB 2 - C22C21 + C11C21 - 2A{D 



(5.36) 
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The eigenvalues of — can then be determined to be : 
T til n + m 33 1 r— — — — — - — 

\{ 2 = + - v(m 11 - m as) 2 + 4m i 3 m 2 4 - 4m M m 2 s 

= m 11+ 2 m 33 - - vT^ ii - ta 33 ) 2 + 4m 13 m 24 - 4m 14 m 23 

(5.37) 



The logarithmic negativity then becomes 

C M {p) = -2 (log 2 (min(l, |A£ 2 |)) + log 2 (min(l, |A£ 4 |))) (5.38) 



Figure 5.1 shows the time evolution of the logarithmic negativity for three 
values of s. The sharp decrease to zero is clearly visible, showing, as one would 
expect, the disentanglement between the particles as they are placed in their 



respective baths. Figure 5.1 also shows that for a constant d, the greater the 
s, the faster the loss of entanglement. 

One may interpret our initial state as a correlated pair of wavepackets 
with width d and the distance between them s. The entanglement then can 
be understood as the interference between the packets. At a distance s = 2d, 
the interferences are destructive and thus there is no entanglement. At any 
other distance, the interference pattern is more or less well-defined which is 
described by a certain value of entanglement. This picture is rather crude, 
yet serves quite well to illustrate why the entanglement would be lost more 



quickly when s increases. If one recalls Figure |2.1| where the entanglement 
entropy is plotted as a function of s, one may recall that as s = 2d, the entan- 
glement disappears, but that it is present for smaller and for larger s, though 
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in smaller amount as s increases. On Figure 5.2 one can also observe that 
as the time increases, the range of s around 2d at which the entanglement 
vanishes increases. This suggests that as time increases, the wave-packets 
would already spread so that the distance around 2d at which the interfer- 
ences become destructive is "blurred". It also suggests that entanglement 
will be present at large s, though in lesser amount still as time increases. 



Dodd and Halliwell DH04 studied disentanglement arising in a separated 
system and an EPR pair in a similar settings but with negligible dissipation 
and using merely a separability criterion, whereas our study concerns entan- 
glement itself. 
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Figure 5.1: Logarithmic negativity vs t for three values of <• 
The values of s are : blue : s = 0.25, red : s = 1, black : s = 
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Logarithmic Negativity for free evolution vs s 



9 



gm1 = 


1,gm2 = 


1,t= 0.001 


gm1 = 


1,gm2 = 


1 , t= 0.005 


gm1 = 


1 , gm2 = 


1 , t= 0.01 




s 



Figure 5.2: Logarithmic negativity vs s for three times 
The plot are for : blue : t = 0.001, red : t = 0.005, black : t = 0.01 



84 



Harmonic Potential Between Two 
Particles in a Heat Bath 



The previous chapter has allowed us to observe that the entanglement 
between the two particles disappears quickly when the particles are free ex- 
cept for their interaction with the heat bath. In the present chapter we will 
now look at the influence of a harmonic potential between the particles. We 
show that allowing the particles to interact may delay the vanishing of the 
entanglement. We also find a striking difference in behaviour between the 
over and the under-damped cases. 



6.1 Time Evolution with a Harmonic Potential 

We would like to see how allowing the two particles to interact may 
influence the evolution of the entanglement and maybe slow down its de- 
crease. For that, we add a harmonic interaction between the two-particles. 
A bipartite harmonic potential Hamiltonian can be written as 



p\ p\ mul 2 
ff .= 2^ + 2^ + ^-(^-^) (6.1) 
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where the r ^ s -{x\ — x 2 ) 2 represents the interaction. 

With position coupling, the master equation becomes 



P = ~\ [H a ,p\ ~jr [xi, [xi,p] + ] - [xi, [x u p]] 



-lh 1 [«]J-f[, 1 |„I (6.2) 
Writing the density matrix in position representation, p(xi, x 2 ; 2/1,2/2), we get 



dp ih ( d 2 d 2 d 2 & 

+ -7^2 )P 



dt 2m \dxl dy\ dx\ dy 2 

{(xi - x 2 ) 2 - (y 1 - y 2 ) 2 ) p 



imui 2 . 



2h 

71 / x ( d d 
2m 

72 /_ .. x ( 9 d \ 7 2 A;T 2 



(d d\ lx kT x 



2m {X2 - y2) \dx- 2 -dyT,) P ~^ {X2 ~ V2) P (6 ' 3) 



If we now perform the customary change of variables x = u + hz,y = u — Kz 
and p(x, y, t) — > P(u, z,t), then 

(xi - x 2 ) 2 - (yi - I/2) 2 —(ui + hzi —u 2 — hz 2 ) 2 - (ui - hzi -u 2 + hz 2 ) 2 

=4h{u 1 -u 2 ){z 1 -z 2 ) (6.4) 

so that 



dP r 1 ( d 2 d 2 \ 

( u > z >*) = {tt-I l^r + 7^r) -2imu 2 (u 1 -u 2 )(z 1 -z 2 ) 



dt ' 1 2m \duidz\ du 2 dz 2 

7i d A 2 72 9 



- — z 1 — - 4 7 iA;Ti ^ 2 - ^ z 2 — - 4 72 A;T 2 2 2 2 )p(u, z, i) 
m ctei m az 2 J 

(6.5) 
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Let us apply the Fourier Transform 



P(q,z,t) = J d Ul du 2 P(u,z,t)e- iqiUl - iq2U2 (6.6) 



whose inverse is 



P(u, z,t) = A y dgi rfg 2 P(q, z, t) e wi«i-H«t« (6 7) 



Since 



- i( Ul - u 2 )Pe- iqiul e~ mU2 = [ A _ A j p e -wi«i e -wa«2 ( 6 g) 

1 9gi 9g 2 / 



we get the differential equation 
dP. 

^(q,z,t) 

I f d 9 \ 71 9 72 9 

'o~~ ?1 a - + 92 7T~ 21 » Z2 J^- 
lm \ ozi oz 2 J m oz\ m oz 2 

+ 2mu 2 ( A _ JL\ ( Zl _ Z2 ) _ 4 7lA;ri ^ _ 472A;T2 ^2 ^ p(q> z t) (6 Q) 



This equation can be solved using the method of characteristics. We have 
the characteristic equation 

<9v M , 

TT- = ^ V (6.10) 

dt 2m y J 



with v T = (z\, z 2 ,qi,q 2 ) and 
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/ 



M 



2 7l 





1 








272 





1 


-4m 2 Uo 


Am 2 ojl 








4m 2 o)Q 


-4m 2 Wo 









On a characteristic we then have 



(6.11) 



dP 
~dt 



(6.12) 



Taking 71 = 72 = 7 and T ± — T 2 — T for simplicity, the eigenvalues and 
eigenvectors of M can be computed to be 

A T = (0, 2 7 , 7 + yV-Sm 2 ^ 2 , 7 - yV-Sm 2 ^) = (Ai, A 2 , A + , A_) 

(6.13) 



and 



Q 



\ 



1 


1 


1 


1 


'27 




A_ 


A+ 


1 


1 


1 


1 


'27 




A_ 


A+ 


1 





— 1 


-1 


1 





1 


1 



(6.14) 



/ 



Since Q 1 MQ = D where D is the diagonal matrix, we need Q 1 : 



Q- 1 









1 

2 


1 

2 


1 


1 


J_ 


_1_ 


2 


2 


47 


47 


A+A_ 


A+A_ 


A_ 


A_ 


2(A+-A_) 


2(A+-A_) 


2(A+-A_) 


2(A+-A_) 


A + A_ 


A + A_ 


A+ 


A+ 


2(A+-A_) 


2(A+-A_) 


2(A+-A_) 


2(A+-A_) 



\ 



(6.15) 
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Using these results, one can rewrite the differential equation as 



2m-^ =QDQ-\ 
at 



dw D 



w 



dt 2m 
with w =Q~\ 



(6.16) 
(6.17) 
(6.18) 



This is easily solved : 



or more explicitly 



We can then write 



with 



Dt/2m 



w(t) = w(0) e Dt/2m 



w 1 (t)=w 1 (0) 
W2 (t) = w 2 (0)e /t/m 
w 3 (t) = w 3 (0)e x +^ 2m 
W4 (t) = W4 (0)e A -*/ 2m 

v(f) = Qe Dt / 2m Q-\ Q 



(l 

e< t/m 

















\ 



(6.19) 



(6.20) 



(6.21) 



e \-t/2m 
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We get 



w = Q *v 



2(A+-A_) 

-A+A_jziQ+A+A-Z2Q— A+gi +A+g2q 
2(A+-A_) 









1 

2 




1 
2 


1 


1 


J_ 




J_ 


2 


2 


47 




47 


A+A_ 


A+A_ 


A_ 




A_ 


2(A+-A_) 


2(A+-A_) 


2(A+- 




2(A+-A_) 


A + A_ 


A + A_ 


A+ 


A+ 


2(A+-A_) 


2(A+-A_) 


2(A+- 


-A-) 


2(A+-A_) 


2 


_l_ 92 

T 2 








ZlO 1 z 20 
2 ~ r 2 


, 910 _i_ 92 

T 47 T 47 








A_|_A_2:iq — A_|_A_ 


Z20+A-91Q - A_ 


92 







\ I z \ 

^10 



^20 

9io 

\ 920 / 



(6.22) 



Then 



/ 





3 7*/ m 






zip 
2 







gA + t/2?Tt Q 

Q e \-t/2m 



9lp I 92 
2 ~ r 2 



I 220 1 910 1 920 

2 47 47 



) 



^10 
2 



91 
2 

220 



+ 



92p 
2 



\ 



I 910 I 92Q 
47 47 



A+ A_ — A + A_22o+A_gi — A-ggp 
2(A+-A_) 

— A+A-Z1Q+A+A-Z2Q — A+gip+A-f 320 
2(A+-A_) 



3 7*/m 



j e A+t/2m 



A + A_zi -A + A-Z2o+A-(?i -A-g2o 
2(A+-A_) 

-A + A-zi +A + A-Z2o-A + gi +A + <72 \ \_t/2m 

2(a+-a_) ; e 



(6.23) 



and 



v (t) =Q e Dt / 2m w 
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1 


1 


1 


1 


'27 




A_ 


A+ 


1 


1 


1 


1 


'27 




A_ 


A+ 


1 





-1 


-1 


1 





1 


1 



We finally get v as 



glQ I 92 Q 
2 ' 2 



( £io i £20 _i_ 9io _i_ £20 1 p 7*/ m 
^ 2 ^ 2 T 4 7 T 4 7 j C 



j e A + t/2m 



A + A-Zi -A + A-Z2Q+A-gi -A-<3'2() 
2(A+-A_) 

^A + A-Zio+A + A-22o-A + qi +A + g2 \ A_t/2m 

2(A+-A_) ; e 









Z2 






Qi 










V ^2 / 



where the matrix J-(t) is given by 



(6.24) 



( a+(t) a-(t) 8+(t) 8-(t) \ 

a-(t) a+(t) <L(t) 6+(t) 

-Pit) Pit) u + (t) u_(t) 

\ p(t) -P(t) i/_(f) u + (t) J 



with 



a±(t) =— ± ! 



(6.25) 



Pit) =A + A_ 



2 " 2(A+ - A ) 

e A + t/2m _ e A_t/2m 



ai(t) ± a 2 (t) 



2(A+-A_) 



■y g7</m gA+t/2m gA_t/2m 

5±(t)= -4^ + ^ ± 2(A + -A_) 



5i(t)±5 2 (t) 



1 A e A_t/2m _ A _ e A +t /2r, 

^=2* 2(A + -A_) 



1 



±1/(0 



(6.26) 
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We can then insert zf(t) and z\{t) into (6.12), we get 



^ = -A 7 kT(z! + 4)P 

= -4 7 A;t| [zxl + z 2 2 ) (a + {tf + a.{tf) + 4z 10 z 20 a + (t)a_(t) 

+ (9io + ?2o) (5+(t) 2 + 5-(t) 2 ) +4q l0 q 20 S4t)5 + (t) 
+ 2f>i gio + z 20 q 2o ) (a_(f)<L(t) + a + (t)5+(t)) 

+ 2(z 10 q 20 + z 20 q 10 ) («_(*)$+ (t) + f> P ( 6 - 27 ) 

After simplification, we can write 
dP 

Tt = - 4 ^ T< 

, o ,Je^l m \ 2 + e x+t / m \2 e x.t/ m \ + \_e<t/r 



+ 4:Z 1Q Z 20 



2 2(A+-A_) 2 2(A+-A_) 2 (A+ - A_) 

,2-yt/m )2 e X+t/m ^2 gA_i/m A i A_e 7 */ r 



4 4(A+-A_) 2 4(A+-A_) 2 2(A+ - A_) 2 
+ (gio + Q2I) 

I g27t/m g7t/m gA+t/m gA_t/m ^yt/m 

X I — + — -= + — ~ T77 + 



^2 g 7 2 4 7 2 2(A+-A_) 2 2(A+-A_) 2 (A+ - A_) : 
+ 4gi g 20 

]_ e 27i/m e 7*/ m gA+t/ m gA_t/m g7*/ m 
y I I l_ 



I67 2 I67 2 8 7 2 4(A+-A_) 2 4(A+-A_) 2 2(A+ - A_) 2 

+ (*10910 + 2 2092 ) 

( e 2-yt/m e 7*/ m A + e A+ * // " 1 A ^— t l m 2^e lt ^ m 
2^ + (A+ - A_) 2 + (A + - A_) 2 " (A + - A_) 
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+ (^10<?2 + ^209l ) 

( p 2~ft/m p-yt/m \ p \+t/m \ X-t/m 9 yt/m 



This can easily be integrated : 

P = Pq exp [-4rykT (xi(zi 2 + 220) + #1^10^20 + X2(qil + 420) + ^2?i ?2 )] 
x exp [-47/cT (Ai(zio?l + ^20<?2 ) + A 2 (2i g2 + ^20?lo))] 

with 

Xl 4 7 l j + 2(A + -A_) 21 j + 2(A+-A_) 21 j 



7 (A + -A_)^ 

1 2 7 l j ( A+ -A_)2 l J (A + -A_)2 le ij 

16m^_ _ 
7(A + -A_)^ e iJ 

X2 8 7 2 16 7 sl J 4 7 3 le j + 2(A + -A_) 2 A + 

m ( e A_t/m _ ^ m ( e 7t/n» _ 1) 

+ 2(A + -A_) 2 A_ " 7 (A + -A_) 2 

t + iIL( e 27*/m_ 1 )_iIL( e 7t/m_ 1) m(e A +*/ m - 1) 



4 7 2 8 7 3V 7 2 7 3V ' (A+-A_) 2 A H 

_ m(e A - t/m - 1) 2m(e /t/m - 1) 
" (A + -A_) 2 A_ + 7(A+-A_) 2 



Ai = — ( e 2 ^ m -1)- — (e^ t/m -1) + 
1 4 7 2 le ij 2 7 2 le ij+ (A + -A_) 2 ' (A+-A_) 



' 2 7 2 ^ 
2m(e^ m - 1) 
(A + -A_) 2 

A 2 = i!L( e 27t/m _ 1} _ _^_( e 7*/m _ x) ^(e A+ * /m - 1) m(e A -*/ m - 1) 



2 



4 7 2V ' 2 7 2V ' (A + -A_) 2 (A + -A_) 

2m(e^/™ - 1) , 

+ (A + -A-)' (6 ' 29) 

93 



We recall the initial state (5.14) 



P(q , z ; 0) = exp [-e + h 2 Zil - e +^2o + 2e -^ ^10^20] 
x exp 



4(e 2 - e 2 



4(4 iei) 9l o 2 (el-el) qi ° q2 ° 



Using v = -7 r (— t)v, we can readily get 



( z \ 




( Z \ 




= H-t) 


Z2 


QlQ 




Qi 


\ <?2 j 




U 2 / 



Then we can write 



(6.30) 



P = exp [-(e+fr 2 + 4 7 fcT X i) (a+(-t) Zl + a + {-t)z 2 + 6+{-t) qi + 5-(-t)q 2 ) 2 ] 
x exp [-(e+ti 2 + A-fkTxi) (a + (-t) Zl + a+(-t)z 2 + <J_(-t)?i + M - *)?*) 2 ] 
x exp[(2e_fr 2 - AjkTOt) 

x (a+(-i)^i + aq_(-i)z 2 + + M _ *)92) 

x (a+(-t)zi + a + (-t)«2 + $-{-t)(-t)qi + &+.(— *)g 2 )] 



x exp 
x exp 



'4(e 



-( 



4(4 - e 2 



+ A 1 kT X 2) (~P{-t)zi + P(-t)z 2 + v+(-t) qi + y-{-t)q 2 y 



+ 4 7 A;Tx 2 ) (P(-t)zi - 0{-t)z 2 + v-(-t)q x + v + {-t)q 2 )' A 
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xex PH 2(e 2 + -6 2 ) +^kT9 2 ) 

X (-p(-t) Zl + P(-t)z 2 + v+(~t)qi + V-(-t)q 2 ) 
x (p(-t) Zl - P(-t)z2 + M-*)?i + "+(-*)92)] 
x exp[-47/cTAi (a + (-t)2i + a+(-t)22 + M~*)9i + M - *)**) 

X + P(-t)z2 + + M"*)**)] 

x exp[-47/cTAi (a+(-t)zi + a+(-t)22 + + 5+(-t)g 2 ) 

x (p(-t)z 1 - f3(-t)z 2 + + v+(-t)Q2)] 

x exp[-47/cTA 2 (a + (-t)z! + a + (-t)z 2 + <J+(-t)9i + <5— ( — 
x (/9(-t)*i - /0(-t)*2 + M-*)?i + f+H)**)] 

x exp[-47/cTA 2 (a + (-t)z! + a + (-t)z 2 + <L(-t)gi + <5-h( — 

x (-P(-t) Zl + P(-t)z2 + + M"*)^] (6.31) 

After some more unpleasant algebra, we can write 

P = exp [— Aq\ — Aq\ — £qiq 2 — Bz\ — Bz 2 — T>Z\Z 2 
- C x z x q x - C x z 2 q 2 - C 2 z x q 2 - C 2 z 2 q x ] 

with 

A =(5i(-t) 2 + 5 2 (-tf)(e + h 2 + 4 7 A;T Xl ) - (S 1 (-t) 2 - 5 2 (-tf)(2^h 2 - A^TO,) 

+ (\ + ,/ H) 2 )( 4 (4_ e 2 ) + W») + (| - y H) 2 )( 2 (4- e 2) + 4 ^ 2 ) 

+ 4 7 A;T ((Sii-t) + 26 2 (-t)is(-t))A 1 + (^(-i) - 25 2 {-t)v{-t))\ 2 ) 

(6.32) 
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(«i(-t) 2 + a 2 {-tf ){e + h 2 + 4 7 A;T X i) - {a^-tf - a 2 {-t) 2 )(2e_h 2 - 4 7 /cT#i) 
+ 2/3(-t) 2 ( 4( ^+ £2) + 4 7 A;T X2 ) - /3(-t) 2 (^_^ + 4 7 CT 2 ) 

+ &ykTp(-t)a 2 (-t) (A 2 - Ai) (6.33) 

4(a!(-t) 2 - a 2 (-t) 2 )(e + ft 2 + 4 7 A;T X i) - («i(-t) 2 + a 2 (-t) 2 )(2e_ft 2 - 4 7 A;T^) 
" 4/3(-t) 2 ( 4( ^+ £2) + 4 7 A;T X2 ) + 2/3(-t) 2 (^— ^ + 4 7 A;Tfl 2 ) 

+ l&ykTP(-t)a 2 (-t) (Ai - A 2 ) (6.34) 

4(5!(-0 2 - 5 2 (-t) 2 )(e+n 2 + 4 7 A;T X i) - (5i(-t) 2 + 5 2 (-tf)(2e_h 2 - 4 7 A;T^) 
+ 4 (| " K-*) 2 )( 4(e2 ^ e2) + 4 7 A:T X2 ) + (i + K-t) 2 )( 2(g2+ e : e2) + 4 7 A;m) 

+ 8 7 A;T - 2(J 2 (-t)^(-0)Ai + (Si(-t) + 2tf 2 (-i)i/(-t))A 2 ) 

(6.35) 

4(a 1 (-t)5 1 (-t) + tt2 (-t)5 2 (-t))( e+ n 2 + 4 7 A;T X i) 
- 2(a 1 (-t)5 1 (-t) - a 2 {-t)5 2 {-t))(2t^h 2 - kykTOx) 

+ 4/3(-t)»/(-t)( 4 C + g2) + 4 7 A;T X2 ) - 2/3(-t)u(-t)( ^=- ^ + 4 7 A;T£ 2 ) 

+ 4 7 A;T( ai (-t) + 2a 2 (-t)z/(-t) - 2/3(-*)<5 2 (-t))A 1 

+ 47fcT(ai(-t) - 2a 2 (-t)i/(-t) + 2/3(-t)5 2 (-t))A 2 (6.36) 

4(ai(-t)(Ji(-i) - a 2 {-t)5 2 {-t)){e + h 2 + 4 7 A;T X i) 
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- 2(a 1 (-t)S 1 (-t) + a 2 (-t)5 2 (-t))(2e_h 2 - A^kTOx) 

- Af3(-t)v(-t)( + 4 7 fcT X2 ) + 2/3(-t)K-t)( 2 £ ~ + 4 7 m? 2 ) 

+ 4 7 fcT(« 1 (-t) - 2« 2 (-t)z/(-t) + 2f3(-t)5 2 (-t))A 1 

+ ^kT( ai (-t) + 2a 2 (-t)z/(-t) - 2f3(-t)5 2 (-t))A 2 (6.37) 



Covariance Matrix and Logarithmic Negativity 



Using (5.17) and (5.28 - 5.31), the covariance matrix can then be written as 



3 



AA -Ci 2S -C 2 

-d B -C 2 V/2 

28 -C 2 AA -Ci 

-C 2 V/2 -d B 



(6.38) 



The partial transpose is then 



3 



AA 




28 


-c 2 


C x 


B 


c 2 


-V/2 


28 


c 2 


AA 


-Cx 


-c 2 


-V/2 


-Cx 


B 



(6.39) 
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One can again calculate — 



1 T 1 



to get 



T — T 



11 



( tl 11 tl 12 tl 13 tl 14 ^ 

tl 21 tl 22 tl 23 tl 24 

tl 31 tl 32 tl 33 tl 34 

y tl 41 tl 42 tl 43 tl 44 J 



where tl 12 = tl 21 = tl 34 = tl 43 = 



tl 11 = tl 22 = tl 33 = tl 44 =4AB -VS + C 2 2 - C\ 

tl 13 = tl 24 = tl 31 = tl 42 =2SB - 2AV 

tl 14 = -tl 32 =CiV - 2C 2 B 

tl 23 = -tl 41 =8AC 2 - 4ECi 



The eigenvalues of — T °"~^ T can tnen ^ e determined to be 



Xl >2 =tl 11 + ^/tl 13 - tl 14 tl 



23 



AL =tl 11 - Jtl ? 3 - tl i 4 tl 2 :; 



The logarithmic negativity then becomes 



C N {p) = -2 (log 2 (min(l, |A£ 2 |)) + log 2 (min(l, |Aj 4 |))) 
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6.2 Observations and Remarks 



Figure 6T shows how the harmonic potential influences the entangle- 
ment. One can notice that in the highly over-damped case (7 = 3 and 
Uq = 1), the entanglement vanishes as quickly with the potential as it does 



without the potential. Figure [672] illustrates the over-damped behaviour. One 
may easily notice that all the curves seem to coincide, suggesting that an 
over-damped harmonic interaction may do little to improve on ESD. More- 



over Figure 6J3 shows that in fact, the entanglement vanishes at shorter 
times, though exponentially close to the free evolution vanishing time. In 
the slightly under-damped case (7 = 1.5), the entanglement also disappears 
around the same time as in the free evolution case. In the highly-under- 
damped case (7 = 0.2), one can observe that C decreases non-uniformly, 
disappears then re- appears for a short while. This leads one to wonder how 
keeping the system slightly or more under-damped may help the entangle- 



ment. Indeed, Figure 6.4 lets us observe that as the damping decreases, the 



logarithmic negativity vanishes at longer times. However, if the damping 



decreases further, as can be seen on Figure |6.5[ C oscillates to a constant 



value greater than 0. Figure 6.6 shows more example of this behaviour. 



To understand better this behaviour, one may recall M's eigenvalues, 
namely 



0, 2 7 , 7 + sJ^-Sm^l 7 - sJ^-Mui) = (A 1; A 2 , A+, A_) 
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and notice that the critical eigenvalues A + and A_ become complex as 7 2 < 
8m 2 u 2 . This brings an oscillatory term into the covariance matrix terms and 
consequently, into the symplectic eigenvalues and the logarithmic negativity, 
resulting into the entanglement's oscillatory convergence towards a constant. 

One may recall our physical interpretation of the initial state. The har- 
monic potential may now be interpreted as a pulse, such as that used to 
increase the quality factor in lasers. In the over-damped case, the pulse con- 
flicts with the wavepacket's original interference, resulting in a faster loss of 
entanglement. In the slightly under-damped case, this oscillatory behaviour 
begins to resonate with the wavepacket, sustaining the interference pattern 
for a while. In the highly under-damped case, this resonance dominates and 
the initial interference pattern becomes drowned into a larger interference 
packet resulting from the harmonic potential's oscillations. 



In FT06 , Ficek and Tanas study a two qubits system coupled to a 
radiation field where they allow spontaneous decay of the atoms. They show 
that the entanglement vanishes but is revived twice, with different reasons 
for each revival. The first revival is due to the regaining of coherence due to 
the spontaneous emission, while the second is related to the asymmetric state 



population. In FT08 , the authors study the emergence of entanglement be- 
tween two initially non-entangled qubits due to spontaneous emission, pro- 
vided both atoms are initially excited and in the asymmetric state. They 
show this creation of entanglement to be a function of the separation be- 
tween the atoms. Their results suggest that allowing an interaction between 
two particles initially entangled will delay the vanishing of the entanglement 
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and revive it, or create entanglement between two initially non-entangled 
particles. We show that when using a harmonic potential as the interac- 
tion, the entanglement revival depends on how strong the coupling is with 
respect to the oscillator's frequency. In fact, we show that if the damping is 
sufficiently low, the entanglement survives for very long times. 
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Logarithmic negativity with and without potential 



4 




Figure 6.1: Logarithmic Negativities with and without potential for loq — 1 

The dashed lines represent C plotted without the potential. The plots are 
(7 — 7i — I2)' green : 7 = 3, red : 7 = 1.5 and blue : 7 = 0.2 
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Logarithmic negativity in the highly over-damped case 



0-1 


, o= 0.05 


— 0-1 


, o= 0.1 


0-1 


, o= 0.15 


gm1 


. 1, gm2 = 1 



0.05 - 

- 



0.05 0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 



Figure 6.2: C in the highly over-damped case 

The plots are obtained keeping 7 = 7i = 72 = 1 and letting Uq vary as : 
dashed : Uq = 0, blue : Uq = 0.05, red : Uq = 0.1 and black : Uq = 0.15 
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Logarithmic negativity in the highly over-damped case 







0-1 


, o= 0.05 


4.5 




— 0-1 


, o= 0.1 






0-1 


,0= 0.15 






gm1 


. 1. gm2 = 1 



ra 2 




0.317 0.3172 0.3174 0.3176 0.3178 0.318 0.3182 0.3184 0.3186 0.3188 0.319 



lure 6.3: Logarithmic negativity in the highly over-damped case 



This is a detail of Figure 6.2 and has the same legend 
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Logarithmic negativity in the slightly under-damped case 




Figure 6.4: Logarithmic negativity in the slightly under-damped case 

The plots are again obtained while fixing 7 = 71 = 72 = 1 and letting u; 
vary as : dashed : u = 0, blue : u = 0.2, red : u = 0.5, black : u = 0.7 

and green : 0.9 
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Figure 6.5: £ as the damping decreases 

The plots are obtained with 7 = 1 and : blue : Uq = 0.9, red : ojq = 1, black 

: u = 1.2 and green : ojq = 1.5 
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Figure 6.6: C in the highly under-damped case 

The plots are obtained with 7 = 1 and : blue : u = 1.5, red : Uq = 1.7, 
black : Uq = 2 and green : Uq = 2.2 
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Setting x' = x — fi, we can obtain the off-diagonal variances A/if and 
Afi 2 , and study the decoherence. If we recall x = u + hz and x' = u — hz, we 
can notice fi = 2hz. Then 

(/if) = J /iip(0, x', t)dxi dx 2 
fi 2 P(0, z, t)dzi dz 2 
-Ah 2 J zfP(0,z,t)dz 1 dz 2 
-Ah 2 J z\e- Bz "' Bz -- Vz ^dz x dz 2 



Ah 2 



Alt' 



2B 



AB 2 - V 2 AB 2 - V 2 



(6.43) 



We can easily see that (hi) vanishes as 



(/ii) = J jj 2 p(0,x',t)dxidx 2 

=Ah 2 J Zie -^l-v4-vz, Z2 dzi dz2 



(6.44) 



Thus we get 



Similarly 



lQh 2 ix 



(AB 2 - V 2 ) 



3/2 



(6.45) 



(6.46) 



Figure |6.7| shows how the coherence inside the system evolves alongside 
C in the highly under-damped case. The variance was rescaled by a factor 
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of one third. It can easily be seen that the coherence peaks as the negativity 
dips. This, however, can be seen not to happen in the slightly under-damped 



case, as shown on Figure 6.8 (note that this time, the variance is resized by 
a factor of one tenth). It is important to note that the coherence in this 
case never vanishes, but converges towards its minimal value at roughly the 
same time the logarithmic negativity vanishes. One may recall Chapter 4's 
results and observe that the coherence in both the slightly and the highly 
under-damped case act in a similar fashion to that of the one-particle case. 

It is quite remarkable that as C converges, so does the coherence but with 
opposite oscillations. This suggests that as the entanglement dies, the state 
regains some coherence, then as it regains coherence, the system regains some 
entanglement and the entanglement is high as the lowest coherence. 
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Figure 6.7: C and Coherence in the highly under-damped case 
The dashed line represent the coherence. The variance is rescaled by 1/3. 
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Logarithmic Negativity and Coherence in the slightly under-damped case 



g= 1,0=0.9 

Cce g= 1 , o = 0.9 



0.5 1 1.5 2 2.5 3 

t 



3.5 4 4.5 



Figure 6.8: C and Coherence in the sightly under-damped case 
The dashed line represent the coherence. The variance is rescaled by 1/10. 
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Conclusion 



This work set out to study how the entanglement in a bipartite Gaus- 
sian state evolves with time, when the state is subjected to an environment. 
To this end, a master equation approach was chosen and the entanglement 
studied. Thus in Chapter 2, the Von Neumann entropy was obtained for an 
initial bipartite Gaussian state. It was then formally shown to be invariant 
under closed system dynamics. A method of studying Gaussian state was 
introduced, namely covariance matrices, which make it very easy to estimate 
the entanglement since it is directly related to the density matrix. Then the 
logarithmic negativity can easily be obtained using the symplectic eigenval- 
ues of the covariance matrix. In Chapter 3, a Non-Rotating Wave master 
equation was derived for a general system Hamiltonian using the Quantum 



Langevin Equation as derived in L088 and a perturbation technique. The 
master equation obtained is similar to that of Savage and Walls. 

In Chapter 4, a single particle state was evolved using the N.R.W. master 
equation in the case of two different Hamiltonians. Results obtained with a 
free-particle Hamiltonian in the first section and with a hamonic potential 
in the second section allowed us to verify the results obtained by Savage and 



Walls SW85b, SW85a . Studying the off-diagonal elements of the resulting 
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density matrix, we observed two different damping behaviours, depending on 
the coupling constant 7 and the frequency u of the oscillator. 

Still using the N.R.W. master equation, the subsequent chapters were 
dedicated to the study of the entanglement in a two-particle Gaussian state. 

In Chapter 5, the bipartite system was coupled to two independent heat 
baths, one for each particle. The independence of the baths ensured that no 
entanglement would be created from the interaction between the particles 
and the baths. This resembles a crudely simple set-up where information 
may be coded in a pair of correlated wave-packets that is afterwards shared 
between two independent parties. 

A free-particle evolution revealed that the entanglement vanishes at very 
short times. In Chapter 6, a harmonic potential was added between the 
two particles to obtain some quite fascinating results. If the systems are 
over-damped, i.e. if the coupling is strong compared to the frequency of 
the oscillator, the entanglement behaves similarly as it would without the 
potential, vanishing at exponentially shorter times. If the systems are slightly 
under-damped, the entanglement decreases at longer times. In the highly 
under-damped case, the logarithmic negativity converges towards a constant 
in an oscillatory manner. 

This result suggests that under certain conditions, it is possible to main- 
tain entanglement in a system for a long time. One can them imagine a 
set-up where two parties share an entangled state in which some informa- 
tion has been encoded. Then the information contained in the entanglement 
would endure. However, the study concerns one type of environment, one 
particular inital state and one particular interaction. It may be that other 
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types of baths will prove less dissipative or that other types of interactions 
or initial states will prove more useful in practical applications. 
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Appendices 



A.l Extra derivations 



Thermal equilibrium relations 



The relations (qjqk), (PjPk), (QjPk) will here derived. The bath Hamil- 
tonian is 



Pj i 2 2 

in , ' ' 



(A.l) 



It can be written in terms of the creation and annihilation operators. 



aj = 



mjUJjqj + ip s 



t _ rrijUjqj - ipj 
a j ~ 



(A.2) 



so that 



h 



Qj = 



2rrijU!j 



[aj + at) 



Pj = 1 



hcOiin 



J J {a)-a 3 ) (A.3) 
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Then 



Pi o 9 HujjTYlj , f .9 n f\j 

rrij 33 2m j 3 3 2oj j m j 



a j + 2a j a) j + a l j • r/j • 2<t a t - r/ ; 



aj + djj 

f 21 



ha 



=hu>j(ajOj + a]aj) = hiOj(2a^aj + 1) = hujj(2hj + 1) (A. 4) 



So 



3^ + -) 



(A.5) 



Now 



h 



h 



h Tr 


(2n + l) e -^ fiw i( ft +3) 


2mjUj 


Tr |e-/^j(™+5) 





ft ^ (2n + l) e -^ n e-^ + ^ 



2m i Uj ^— ' e -pfiwjn e -Ph(uj+±) 

J J n=0 



h 



E 



(2n + l) e -^ n 



J „=o 



2rrijix!j 



{2hj + 1) 



(A.6) 



The sums can be express as follows 



E' 



-fitkUjn 



1 



1 _ g-/3W, 



(A.7) 



and 



n 



1 d 



e -f3hujj 



(1 _ e /3^)2 



(A.8) 
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Hence 



{Wk) 



h 



E 



■> J n=0 



h 



+ 



g—fibwjn 

2 ePnuj 



2m j u j \ 1 - e-P^j 



+ 1)5 



h 1 



e -f>i*> s 



Ivrijbjj 1 



coth | ) 



/i 



coth 



V 2 

fa*;,- 



2A;T 



3j* 



(A.9) 



Similarly 



hrrijUJj t 2 hrrijiOj . + +. 

= -—Z- L (-M j -l)6 jk 



hrrijUj ^ r 


-(2n + l) e _ ^'( ft+ 5) 


2 


Tr 







Ojk 



E 

n=0 



-/3^j?i e -/3ft(wj + |) 



jk 



hrrijUJj A (2n + l) e -^ n 

: 9 e -/3^n "?' fe 



n=0 



hrrijUJj 


coth 1 




2 


V 2 


hrrijUJj 


coth | 




2 


v IkT ) 



(A.10) 
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and 



(QjPk) =i\j h ^l^J _A_(( aj + a))(a) + a,j)) 5 jk = t^(a ja ] - a)a.j) 5 jk 

h 

(qjPk) =t-5 jk (A.ii) 



and finally 



(PjQk) =*y h ^r l \j ^Jt((°] + a i)( a i + 4)) 6 i k = l \^ a \ a i - a i a ]) 6 Jk 

h 

(PjQk) =~ l 2 6 i k = -(QjPk) ( A - 12 ) 
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